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ABSTRACT 

Active  flow  control  of  jets  with  Localized  Arc  Filament  Plasma  Actuators  (LAPP  As)  is  conducted  over  a 
wide  range  of  the  fully  expanded  jet  Mach  numbers  (Mj  or  simply  jet  Mach  number).  The  jet  Mach 
numbers  covered  in  the  present  research  are  0.9  (with  a  converging  nozzle),  1.2  (overexpanded),  1.3 
(perfectly  expanded),  and  1.4  (underexpanded)  with  a  design  Mach  number  1.3.  Additionally,  limited 
experiments  are  carried  out  for  an  Mj  =  1.65  perfectly-expanded  jet.  The  exit  diameter  is  2.54  cm  (1  inch) 
for  all  cases  and  eight  LAPP  As  are  equally  distributed  on  the  perimeter  of  a  boron  nitride  nozzle 
extension.  The  jet  spreading  is  strongly  dependent  on  duty  cycle,  forcing  frequency,  and  azimuthal  modes. 
The  performance  of  LAPP  As  for  jet  spreading  is  investigated  using  two-dimensional  particle  image 
velocimetry  (PIV).  There  is  an  optimum  duty  cycle,  producing  maximum  jet  spreading,  for  each  forcing 
frequency.  A  relationship  between  the  optimum  duty  cycle  and  forcing  frequency  is  determined  from  the 
extensive  results  in  the  Mj  0.9,  and  this  relation  is  used  for  all  experiments.  The  effect  of  forcing 
frequency  is  investigated  for  a  wide  range  of  forcing  Strouhal  numbers  (Stop  =  fpD/Ue,  where  fp,  D,  and 
Ue,  are  forcing  frequency,  nozzle  exit  diameter,  and  jet  exit  velocity  respectively),  ranging  from  0.09  to 
3.0.  The  azimuthal  modes  (m)  investigated  are  m  =  0  -  3,  ±1,  ±2,  and  ±4  -  this  comprises  all  modes 
available  with  eight  actuators.  The  performance  of  LAPP  As  does  also  strongly  depend  on  the  stagnation 
temperature  of  the  jet  and  Mj.  The  effects  of  stagnation  temperature  are  investigated  for  1.0,  1.4,  and  2.0 
times  the  ambient  temperature  in  Mj  0.9  jet  for  very  limited  azimuthal  modes  and  Stop.  In  an  Mj  1.65 
perfectly-expanded  jet,  the  control  authority  of  LAFPAs  is  investigated  for  only  m  =  ±1  and  Stop  ^  0.3. 

The  jet  spreading  increases  with  decreasing  duty  cycle  until  the  limit  of  incomplete  breakdown  is  reached. 
Thus,  the  optimum  duty  cycle  is  the  lowest  value,  at  a  given  forcing  frequency,  which  ensures  complete 
breakdown.  Extensive  experiments  in  Mj  0.9  and  1.3  perfectly-expanded  jets  show  that  the  jet  spreading 
is  greatest  at  about  Stop  =  0.3  for  most  azimuthal  modes.  The  most  and  least  effective  azimuthal  modes 
for  mixing  enhancement  are  m  =  ±1  and  ±4,  respectively.  The  results  also  show  that  the  effect  of  forcing 
is  very  similar  in  Mj  0.9  subsonic  and  1.3  perfectly-expanded  supersonic  jets.  The  results  in  the  heated  Mj 
0.9  jet  show  that  the  effects  of  forcing  increase  with  increasing  stagnation  temperature.  In  addition,  the  jet 
spreading  in  m  =  0  is  comparable  to  that  in  m  =  1  at  an  elevated  stagnation  temperature  while  it  was  one 
of  the  less  effective  modes  in  an  unheated  jet.  The  turbulent  kinetic  energy  along  the  jet  centerline  also 
increases  significantly  near  Stop  =  0.3  for  most  azimuthal  modes. 

In  off-design  conditions  of  Mj  =  1 .2  and  1 .4,  the  forcing  is  less  effective  compared  to  the  perfectly- 
expanded  case  of  Mj  =  1.3.  Flow  structure  visualization,  using  Galilean  streamlines,  shows  that  there  are 
naturally-amplified  flow  structures  (generated  by  a  natural  feedback  mechanism  in  the  over-  and  under¬ 
expanded  jets)  in  addition  to  the  structures  generated  by  forcing.  The  competition  of  these  structures 
seems  to  be  responsible  for  the  reduced  effectiveness  of  forcing.  The  performance  of  LAFPAs  in  a  high 
Mj  supersonic  jet  (Mj  =  1.65)  shows  reduced  forcing  effectiveness.  It  has  not  yet  been  conclusively 
determined  if  the  reduction  in  effectiveness  is  due  to  the  lack  of  LAPPA  control  authority  or  increased 
flow  compressibility. 
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1.  INTRODUCTION 

Many  researchers  have  worked  on  jet  flow  control  to  enhance  mixing  and/or  reduce  noise.  Most  of 
the  earlier  jet  flow  control  was  done  in  low-speed  and  low  Reynolds  number  flows.  In  such  flows, 
acoustic  drivers  were  successfully  used  since  the  flow  momentum  and  associated  flow  characteristic 
frequency  are  low.  However,  the  acoustic  driver  does  not  have  sufficient  bandwidth  and  amplitude  in 
high-speed  and  high  Reynolds  number  flows,  since  characteristic  flow  frequency  and  flow  momentum 
increase  as  the  jet  speed  and  Reynolds  number  rise. 

1.1  Jet  instabilities 

The  most  successful  manipulation  of  the  jet  flow  is  related  to  controlling  the  jet  characteristic 
instabilities.  There  are  two  major  instability  modes  in  a  jet:  the  initial  shear  layer  instability  and  the  jet 
column  instability.  These  modes  are  based  on  two  length  scales  in  a  free  jet:  the  initial  boundary  layer 
momentum  thickness  (6)  at  the  nozzle  exit  and  the  nozzle  exit  diameter  (D)  for  a  circular  nozzle  or  the 
nozzle  exit  height  (h)  for  a  rectangular  nozzle.  The  initial  shear  layer  instability  frequency  is  scaled  with 
the  momentum  thickness  (6)  at  the  nozzle  exit.  The  jet  column  instability  or  the  jet  preferred  mode  is  the 
instability  around  the  end  and  downstream  of  the  potential  core,  and  its  frequency  is  scaled  with  the 
nozzle  exit  diameter  (D)  or  height  (h).  The  corresponding  Strouhal  numbers  are  Ste  (=  f6/Uj)  and  Sto  (= 
fD/Uj)  for  initial  shear  layer  instability  and  jet  column  mode,  respectively.  The  f  and  Uj  are  instability 
wave  frequency  and  the  jet  exit  velocity,  respectively. 

The  shear  layer  of  an  unforced  jet  in  the  vicinity  of  the  nozzle  exit  is  very  thin  so  that  its  behavior 
is  very  similar  to  that  in  a  planar  shear  layer,  since  curvature  effects  are  negligible.  The  mixing  layer  near 
the  exit  of  the  jet  is  referred  as  initial  shear  layer.  In  the  initial  shear  layer,  the  maximum  amplification  of 
disturbances  seems  to  occur  around  the  Strouhal  number  (Ste^fQ/Uj)  of  0.012  in  unforced  jets  [Zaman  and 
Hussain  1981],  while  the  maximum  amplification  rate  of  disturbances  occur  around  Ste^O.Ol?  [Freymuth 
1966,  Michalke  1965]  in  forced  jets.  The  input  excitation  amplitude  required  to  control  this  instability  in 
low-speed  flows  is  very  small  and  linear  instability  analysis  has  been  used  extensively  to  explore  various 
aspects  of  this  instability  [Michalke  1965].  When  the  initial  shear  layer  is  forced,  the  increased 
amplification  rate  leads  to  earlier  saturation  of  amplification  and  breakdown  of  amplified  instability 
waves/vortices  into  smaller  scales  so  that  the  amplification  of  instability  is  smaller  than  that  in 
unperturbed  jets  [Zaman  and  Hussain  1981].  Thus,  turbulence  intensity  in  the  downstream  region  can  be 
reduced  when  the  initial  shear  layer  is  forced  at  Ste^  0.017.  However,  the  growth  of  instability  at  Ste^ 
0.012  leads  to  the  large  scale  structures  in  the  shear  layer  of  the  jet,  which  are  responsible  for  the 
entrainment  of  ambient  air  into  the  jet  and  gross  mixing  with  the  jet  fluid. 

The  maximum  amplification  of  the  jet  column  instability  occurs  over  a  wide  range  of  Sto  from  0.2 
to  0.6  [Crow  and  Champagne  1971;  Gutmark  and  Ho  1983;  Ho  and  Huerre  1984;  Cho  et  al.  1998], 
depending  heavily  upon  the  experimental  facility  as  well  as  on  what  is  measured  and  where  it  is  measured. 
This  is  presumably  due  to  the  variations  in  the  naturally  occurring  disturbances  in  the  facilities.  The  jet 
column  mode  can  be  excited  directly  by  forcing  the  mode  with  high  enough  amplitude  [Cho  et  al.  1998]. 

The  initial  shear  layer  instability  and  the  jet  column  mode  can  be  coupled  when  the  boundary  layer 
at  the  nozzle  exit  is  laminar  [Ho  and  Hsiao  1983].  The  coupling  occurs  through  an  integer  number 
(usually  3  or  4)  of  pairings  of  relatively  small  structures  in  the  initial  shear  layer.  Kibens  [1980]  also 
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observed  a  eoupling  of  these  two  modes  in  a  foreed  jet  with  an  aeoustie  driver.  However,  the  Strouhal 
number  along  the  lip-line  of  the  jet  was  not  stepwise,  but  smoothly  ehanged.  This  suggests  that  the  pairing 
did  not  oeeur  in  an  orderly  manner  so  that  the  eoupling  of  the  two  modes  perhaps  did  not  happen 
[Ginevsky,  et  al.  1974].  Based  on  earlier  results  in  Maeh  0.9  subsonie  [Kim  et  al.  2009a]  and  Maeh  1.3 
supersonie  [Samimy  et  al.  2007b]  jets,  it  seems  that  the  jet  eolumn  mode  is  direetly  foreed  by  the 
LAPP  As. 

In  addition  to  the  two  instability  modes  diseussed  above,  there  is  azimuthal  mode  instability  in  a 
eireular  jet.  The  jet  eolumn  instability  is  unstable  to  azimuthal  or  helieal  modes  [Cohen  &  Wygnanski 
1987].  The  major  faetor  deeiding  the  growth  rate  and  amplitude  of  various  azimuthal  modes  seems  to  be 
the  ratio  of  the  jet  nozzle  exit  diameter  to  the  boundary  layer  momentum  thiekness  at  the  nozzle  exit 
(D/6o).  Linear  stability  analysis  of  Miehalke  [1977]  and  Plasehko  [1979]  and  experimental  work  of  Cohen 
and  Wygnanski  [1987],  Corke,  Shakib,  and  Nagib  [1991],  and  Corke  and  Kusek  [1993]  showed  that  for 
large  D/Bq  (D/0  o»l),  both  axisymmetrie  (m=0)  and  the  first  spinning  or  helieal  modes  (m=+l  or  -1)  are 
unstable  in  the  initial  jet  shear  layer.  Linear  stability  analysis  of  Cohen  and  Wygnanski  [1987]  also 
showed  that  for  a  very  thin  boundary  layer  (or  very  large  D/6o),  many  azimuthal  modes  are  unstable  in  the 
initial  shear  layer  region.  Linear  stability  analysis  of  Miehalke  [1977]  also  showed  that  further 
downstream  in  the  jet,  where  the  veloeity  profile  is  bell-shaped,  the  jet  ean  only  support  helieal  modes.  It 
has  also  been  reported  that  the  growth  region  of  helieal  modes  move  further  upstream  towards  the  nozzle 
as  the  jet  veloeity  inereases  [Ho  and  Huerre,  1984].  A  more  detailed  diseussion  on  the  instability  related 
to  a  eireular  jet  is  found  in  Samimy  et  al.  [2007b]. 

1 .2  Heated  jet 

Previous  work  on  eontrol  of  heated  jets  has  foeused  on  eharaeterizing  the  ehanges  in  the  flow  affeeted 
by  varying  the  temperature  of  the  jet.  It  has  been  observed  that  jet  total  temperature  is  an  independent 
parameter  governing  the  ehanges  in  turbulent  boundary  layer  integral  eharaeteristies  (e.g.  displaeement 
thiekness,  ete.)  in  addition  to  jet  exit  Reynolds  number  over  the  range  examined  in  those  experiments 
(280,000  -  1,400,000)  [Lepieovsky  1990].  When  the  boundary  layer  of  a  low  Reynolds  number  (-10,000) 
heated  jet  transitions  from  turbulent  to  laminar,  the  mixing  eharaeteristies  ehange  dramatieally 
[Strykowski  and  S.  Russ  1992].  Lepieovsky  [1986]  experimented  with  aeoustieally  exeited  heated  jets 
eoneluding  that:  jet  sensitivity  to  upstream  aeoustie  exeitation  varies  strongly  with  the  jet  operating 
eonditions,  exeitation  threshold  level  inereases  with  inereasing  jet  temperature,  and  jet  preferred  mode 
Strouhal  number  does  not  ehange  signifieantly  under  varying  eonditions.  In  a  eomparison  of  Maeh  0.3 
and  Maeh  0.8  jets,  it  was  eoneluded  that  the  higher  Maeh  number  jet  aehieved  natural  exeitation  due  to 
heating  and  suggested  that  larger  external  foreing  amplitudes  would  be  required  to  observe  any  exeitation 
effeets  [Ahuja  et  al.  1986].  This  study  also  noted  that  the  trend  observed  in  their  experiments  is 
eontradietory  to  theory  [Ahuja  et  al.  1982]  whieh  states  that  exeitation  effeetiveness  should  inerease  with 
temperature.  Turbulent  boundary  layers  are  less  seleetive  (eompared  to  laminar)  about  optimum 
exeitation  frequeney  [Lepieovsky  1989].  Based  on  these  previous  works  LAPP  As  may  have  an  advantage 
over  previous  aetuators  that  ideally  positions  LAFPAs  to  investigate  exeitation  in  heated  jets.  One 
disadvantage  is  that,  at  least  at  the  present  time,  there  is  no  meehanism  for  direetly  eontrolling  the  foreing 
amplitude  of  LAFPAs  so  an  inerease  in  exeitation  threshold  would  only  be  notieed  if  the  threshold 
exeeeded  the  foreing  amplitude  required  to  eontrol  the  jet. 
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When  the  jet  is  heated  while  maintaining  a  eonstant  stagnation  pressure  (eonstant  Maeh  number),  the 
Reynolds  number  {Re  =  pUjD/ji)  deereases.  The  equations  used  to  ealeulate  the  Reynolds  number  are:  the 
ideal  gas  law,  isentropie  eompressible  flow  relations,  and  Sutherland’s  formula  for  viseosity.  Thus,  for  the 
relevant  values  o^ D  =  2.54  em  and  Mj=  0.9  for  these  experiments,  Reynolds  number  has  the  following 
relationship: 


Re(rJ  =  1.515x10*°  x7;“Vl.302xl0^xr/  (1) 

where  To  is  the  stagnation  temperature  in  Kelvin.  For  example:  Re{290)  -  630,000  and  7?e(580)  -  270,000 
eorrespond  to  a  temperature  ratio  of  about  1.0  and  2.0,  respeetively.  Consequently,  as  the  stagnation 
temperature  rises,  the  Reynolds  number  eould  beeome  suffieiently  small  for  the  boundary  layer  at  the 
nozzle  exit  to  potentially  beeome  transitional  or  laminar. 

Another  way  in  whieh  heated  jets  may  differ  from  unheated  jets  is  the  effeet  of  density  gradients.  The 
density  ratio  between  the  eore  and  the  ambient  air  eould  affeet  the  nature  of  jet  instabilities.  Studies  of 
density  gradient  phenomena  to  this  point  have  foeused  on  mueh  lower  Reynolds  numbers  (typieally  a  few 
thousand)  and  lower  speeds  than  those  in  this  paper,  but  the  eoneepts  may  still  be  applieable.  In  a  jet  with 
suffieiently  large  density  ratio,  there  eould  be  a  region  of  absolute  instability  leading  to  jet  global 
instability.  This  instability  exists  in  addition  to  the  initial  shear  layer  and  jet  eolumn  instabilities  diseussed 
previously.  If  a  suffieiently  large  region  beeomes  absolutely  unstable,  the  jet  may  beeome  globally 
unstable  displaying  oseillatory  behaviors  [Huerre  and  Monkewitz  1990].  Large,  axially  symmetrie 
oseillations  in  the  potential  eore  region  of  jets  with  density  ratio  {S  =  Pjj  Pamb )  (temperature 

ratio  above  1.39)  with  a  eharaeteristie  Strouhal  number  of  -0.3  have  been  observed  [Huerre  and 
Monkewitz  1990,  Monkewitz  et  al.  1990].  Jendoubi  and  Strykowski  [1994]  used  simulations  to  explore 
how  these  instabilities  ehange  with  Maeh  number,  eo-  and  eounter-flow,  and  shear  layer  thiekness.  It  was 
eoneluded  that  inereasing  Maeh  number  deereases  the  level  of  instability.  This  inerease  in  stability  is 
supported  by  the  observation  that  the  region  of  absolute  instability  shrinks  from  S  <  0.7  for 
ineompressible  to  5"  <  0.2  at  Maeh  0.6.  Soteriou  and  Ghoniem  [1995]  reported  on  numerieal  simulations 
of  density  ratios  in  shear  layers  on  the  range  0.33  -  4.0,  whieh  showed  that  as  the  density  ratio  deereases: 
eonveetive  veloeity  slows,  entrainment  inereases,  and  growth  rate  inereases.  In  short,  the  flow  beeomes 
biased  towards  the  denser  fluid.  A  fairly  exhaustive  list  of  additional  studies  ean  be  found  in  the  work  of 
Lesshafft  et  al  [2007]. 

1.3  Objectives  of  the  research 

The  objeetives  of  the  present  researeh  are  multifaeeted  and  related  to  these  questions: 

1)  Does  the  plasma  aetuator,  to  be  deseribed  later,  have  eontrol  authority  in  high-speed  and 
high-Reynolds  number  jets? 

2)  What  are  the  major  parameters  whieh  signifieantly  affeet  the  jet  flow? 

3)  Do  the  aetuators  have  the  same  eontrol  authority  in  off-design  eonditions? 

4)  Are  the  aetuators  more  or  less  effeetive  in  heated  jets? 

5)  What  is  the  range  of  jet  Maeh  numbers  where  the  aetuator  has  eontrol  authority? 

Extensive  experiments  are  earned  out  to  find  answers  for  1)  and  2)  in  Mj  (fully  expanded  jet  Maeh 
number)  0.9  subsonie  and  Mj  1.3  perfeetly-expanded  supersonie  jets  over  a  wide  range  of  foreing 
frequeneies  and  all  available  azimuthal  modes.  To  find  answers  for  3),  experiments  are  eondueted  in  off- 
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design  conditions  of  Mj  1.2  (overexpanded)  and  1.4  (underexpanded).  Based  on  the  results  and 
conclusions  of  Mj  0.9  and  1.3  jets,  very  limited  experiments  are  carried  out  to  answer  questions  4)  and  5) 
in  heated  Mj  0.9  and  unheated  Mj  1.65  jets,  respectively. 
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2.  EXPERIMETAL  FACILITY  AND  TECHNIQUES 


2.1  Nozzle  and  Air  Supply  System 

All  the  experiments  are  eondueted  at  the  Gas  Dynamies  and  Turbulenee  Laboratory  at  the  Ohio 
State  University.  The  eompressed  air,  whieh  is  filtered  and  dried,  is  stored  in  two  eylindrieal  tanks  with  a 
eapaeity  of  43  up  to  16  MPa.  The  eompressed  air  by  three  five-stage  eompressors  is  supplied  to  the 
storage  tank  and  then  to  the  stagnation  ehamber  of  the  jet.  The  air  is  then  diseharged  through  a  nozzle 
with  1.0”  (2.54  em)  exit  diameter.  Three  nozzles  are  used  to  eover  Mj’s  from  0.9  to  1.65.  A  eonvergent 
nozzle  is  used  for  Maeh  0.9  jets,  and  a  design  Maeh  (M^)  1.3  eonverging-diverging  nozzle  is  used  for 
supersonie  jets  in  perfeetly  expanded  (Mj  =1.3)  and  imperfeetly  expanded  eonditions  (Mj  =1.2  and  1.4). 
Shown  in  Fig.  2.1a  is  the  streamwise  eross-seetion  of  1.3  nozzle,  whieh  has  very  a  smooth  eonverging 
seetion  designed  by  the  method  of  eharaeteristies  to  obtain  shoek-free  uniform  exit  veloeity.  Maeh  0.9 
nozzle  has  a  relatively  rapid  eonverging  seetion  eompared  to  the  Maeh  1.3  nozzle  (drawing  is  not  shown). 
For  the  Mj  =1.65  jets,  either  eontoured  or  eonieal  (military  applieation)  nozzle  is  used  as  shown  in  Fig. 
2.1. 


(a)  Md  1.3  contoured  nozzle 


01,00 

i 


(c)  Md  1.65  contoured  nozzle 


(b)  Boron  nitride  nozzle  extension 


01.00 

i 


(d)  Md  1.65  conical  nozzle 


Fig.  2.1  Schematic  of  1.3  and  1.65  nozzles  and  nozzle  extension  (units  are  in  inch). 


At  the  end  of  the  nozzle,  a  boron  nitride  nozzle  extension  is  attaehed  to  house  eight  plasma 
aetuators,  uniformly  distributed  in  azimuthal  direetion  (Fig.  2.1b).  Eaeh  aetuator  is  eomposed  of  two 
tungsten  pin  eleetrodes  with  a  diameter  of  1  mm.  The  eenter-to-eenter  distanee  of  two  eleetrodes  is  about 
4  mm  at  the  tip.  All  eleetrodes  are  plaeed  1  mm  upstream  of  the  extension  exit  within  a  ring  groove, 
measuring  1  mm  wide  and  0.5  mm  deep,  to  prevent  the  plasma  from  being  blown  off  As  shown  in  Fig. 
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2.1,  the  electrodes  are  installed  radially  and  the  tip  of  each  electrode  is  flush-mounded  to  the  inner  surface 
of  the  nozzle  extension. 


Table  1  shows  the  fully  expanded  jet  Mach  number  (Mj,  Mach  number  based  on  the  ratio  of  the 
stagnation  pressure  to  ambient  pressure)  and  corresponding  nozzle  pressure  ratios  relative  to  the  ambient 
pressure.  For  the  M^  1.3  nozzle,  the  jet  is  either  perfectly-expanded  (Mj  =1.3),  over-  (Mj  =1.2),  or  under- 
(Mj  =1.4)  expanded.  The  M^  1.65  nozzles  are  operated  only  at  the  design  Mach  number.  However,  the 
flow  at  the  nozzle  exit  for  the  conical  nozzle  is  not  perfectly  expanded  due  to  its  conical  diverging  section. 
The  effects  of  the  stagnation  temperature  of  the  jet  are  investigated  only  at  Mj  =  0.9  with  stagnation 
temperature  ratios  (To/Ta,  To  =  stagnation  temperature,  Ta=  ambient  temperature)  of  1.0,  1.4  and  2.0.  The 
Reynolds  number  at  the  nozzle  exit  is  also  shown  in  Table  1.  The  Reynolds  number  (Rcd  =  pUeD/p)  is 
based  on  the  nozzle  exit  diameter  (D)  and  velocity  (Ue),  and  also  on  density  (p)  and  dynamic  viscosity  (p) 
at  the  nozzle  exit. 


Table  1  Fully  expanded  jet  Mach  numbers  and  corresponding  nozzle  pressure  ratios. 


Nozzle 

Fully  expanded  jet 
Mach  number  (Mj) 

Nozzle  Pressure  Ratio 
(stagnation  pressure  / 
ambient  pressure) 

Stagnation 
temperature 
ratio  (To/Ta) 

Reynolds 
number  (Rcd) 

Convergent  nozzle 

0.9 

1.69 

1.0 

0.63  X  10" 

1.4 

0.41  X  lO*^ 

2.0 

0.27  X  lO*^ 

Design  Mach  1.3 

1.2 

2.42  (Over-expanded) 

1.0 

1.1  X  10^ 

Converging- 

1.3 

2.77  (Perfectly  expanded) 

1.0 

1.2  X  10® 

diverging  nozzle 

1.4 

3.18  (Under-expanded) 

1.0 

1.4  X  10® 

Design  Mach  number 

1.65 

4.58  (perfectly  expanded) 

1.0 

16.8  X  10® 

1.65 

The  heating  system  is  composed  of  a  Watlow  15  kW  electric  heater  and  a  vertical  heat  storage  tank. 
The  heat  storage  tank  isa3.5m(138  in.)  tall  by  0.3  m  (12  in.)  diameter  cylinder  packed  with  four  sets  of 
vertically  aligned  rows  of  stainless  steel  plates.  An  electric  fan  takes  room  air,  passes  it  through  the 
electric  heat  chamber,  through  the  heat  storage  tank,  and  discharges  it  outdoors.  The  electric  heater  has  a 
maximum  output  temperature  of  866  K  (1 100  T)  which  produces  a  maximum  jet  stagnation  temperature 
of  -775  K  due  to  heat  loss  in  the  storage  system.  During  experiments,  pressurized  air  is  forced  through 
the  heat  storage  tank  to  be  heated  before  entering  the  jet  stagnation  chamber.  The  Mach  0.9  jet 
experiments  can  be  run  continuously  for  approximately  40  minutes  with  minimal  temperature  variation 
-0.2  K/min.  This  system  is  limited,  not  by  a  maximum  flow  rate,  but  by  how  long  the  storage  tank  can 
maintain  a  stable  temperature. 


2.2  Boundary  Layer  Conditions  near  the  Nozzle  Exit 

As  discussed  earlier,  the  condition  of  boundary  layer  near  the  nozzle  exit  plays  a  significant  role 
in  the  jet  development  and  jet  instabilities.  Due  to  the  small  diameter  of  the  nozzle  and  very  thin  boundary 
layer  at  the  nozzle  exit,  even  the  simplest  questions  about  the  boundary  layer  are  nearly  impossible  to 
answer  with  PIV  measurements.  In  order  to  address  this  issue,  a  slightly  larger  converging  nozzle  with 
3.81  cm  (1.5  in.)  diameter  was  used  to  examine  the  boundary  layer  characteristics  of  this  experimental 
setup.  Apart  from  the  change  in  diameter,  the  nozzle  used  in  this  experiment  is  essentially  identical  to  the 
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nozzle  used  for  the  PIV  results.  A  slightly  larger  nozzle  diameter  was  ehosen  to  maximize  the  number  of 
measurement  points  obtained  in  the  boundary  layer  while  remaining  suffieiently  similar  to  the  primary 
nozzle.  Additionally,  a  slightly  larger  diameter  allows  aeeess  to  Reynolds  numbers  (based  on  nozzle  exit 
diameter)  typieal  of  the  operating  eonditions  in  the  unheated  jet  ease  while  avoiding  the  eompressibility 
eomplieations  at  higher  veloeities.  Sinee  hot-wire  measurements  are  very  diffieult  to  obtain  in  a  hot  jet, 
the  deeision  was  made  to  assume  that  any  signifieant  ehanges  in  the  boundary  layer  eharaeteristies  should 
be,  at  least  primarily,  dependent  on  Reynolds  number.  The  free  shear  profile  just  downstream  of  the 
nozzle  was  measured  in  an  unheated  jet  over  a  range  of  Reynolds  numbers  from  200,000  to  600,000 
ereated  by  varying  the  Maeh  number  of  the  jet. 


Fig.  2.2  Normalized  free  shear  layer  velocity  profiles  for  a  range  of  Reynolds  numbers. 


The  normalized  veloeity  profiles  for  the  two  most  informative  eases  are  shown  in  Fig.  2.2.  The 
abseissa  is  the  radial  jet  eoordinate  normalized  by  the  jet  diameter.  Before  proeeeding  any  further,  the 
effeets  of  the  thiek  lipped  nozzle  should  be  noted.  The  eomparatively  large  outer  diameter  of  the  nozzle 
lip  (~3  in.)  ereates  a  reeireulation  region  whieh  aets  to  widen  the  diameter  of  the  jet  as  determined  by  the 
hot-wire  data.  Additionally,  the  wake  profile  typieal  of  thin  lipped  nozzles  is  eompletely  absent.  Without 
a  nozzle  extension  (Fig.  2.2a),  the  normalized  veloeity  profiles  do  not  eollapse  -  indieating  transitional 
behavior  in  the  boundary  layer.  As  Reynolds  number  inereases,  the  profile  pushes  outward  and  the  sloped 
deereases. 

However,  when  a  nozzle  extension  is  attaehed  (Fig.  2.2b),  the  profiles  eollapse  very  well.  The 
only  appreeiable  ehange  with  inereasing  Reynolds  number  is  a  slight  deerease  in  the  eurvature  of  the  high 
veloeity  shoulder.  The  eonsisteney  of  the  profile  is  evidenee  of  a  eonsistently  turbulent  boundary  layer. 
Additionally,  the  extension  slightly  deereases  the  effeetive  jet  diameter.  Through  experimentation  with 
different  nozzle  extensions  (not  shown),  it  was  established  that  the  eollapse  is  eaused  by  either  the  slight 
trip  ereated  by  the  mating  surfaees  between  the  nozzle  and  the  nozzle  extension  or  the  additional  distanee 
provided  by  extension.  In  reality,  both  features  probably  eontribute. 

Sinee  these  profiles  are  of  shear  layers,  not  traditional  boundary  layers,  it  was  deeided  to  estimate 
the  momentum  thiekness  by  fitting  the  profiles  to  a  hyperbolie  tangent  as  performed  by  Beehert  and  Stahl 
[1988].  The  slope  of  the  fitted  profile  is  used  to  as  the  slope  of  a  line.  The  horizontal  distanee  (A)  between 
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where  that  line  erosses  one  and  zero  is  related  to  the  momentum  thiekness  (0)  as  A  =  40.  The  momentum 
thiekness  of  the  LAFPA  extension  ease  was  determined  to  be  -0.09  mm  with  a  variation  of  -0.01  mm 
over  the  range  of  Reynolds  numbers.  The  ease  without  an  extension  had  momentum  thieknesses  ranging 
from  -0.05  to  -0.09  mm.  Previous  work  by  Lepieovsky  [1999]  reports  similar  values  for  the  momentum 
thiekness  of  turbulent  boundary  layers  over  this  range  of  Reynolds  numbers  in  a  5.08  em  (2  in.)  jet.  From 
the  profiles  shown  in  Fig.  2.2b,  the  best  estimate  for  the  boundary  layer  thiekness,  -1.2  mm,  was 
ealeulated  as  the  distanee  from  the  nozzle  lip  {r/D  =  0.5)  to  98%  of  the  free-stream  veloeity. 

2.3  Plasma  Actuator 

The  plasma  generating  system,  shown  in  Fig.  2.3,  has  two  high  voltage  Glassman  DC  power 
supplies,  with  output  of  10  kV  and  1  Ampere.  Eaeh  power  supply  ean  drive  four  aetuators  simultaneously, 
and  thus  up  to  eight  aetuators  ean  be  operated  at  the  same  time.  Eaeh  aetuator  is  eontrolled  independently 
by  a  Behlke  high  voltage  transistor  switeh.  A  National  Instruments  (NI)  analog  board  attaehed  to  a  PC  is 
used  to  generate  eight  independent,  eontinuous  pulse  trains  to  eontrol  the  transistor  switehes.  Details  of 
the  plasma  system  are  provided  in  Utkin  et  al.  [2007]  and  in  Samimy  et  al.  [2007b]. 


< 


Fig.  2.3  Schematic  of  the  in-house  fabricated  8-channel  plasma  generator. 


The  foreing  frequeney,  duty  eyele,  and  azimuthal  mode  are  eontrolled  through  LabView,  NI 
software.  The  available  azimuthal  modes  with  eight  aetuators  are  m  =  0-3,  ±1,  ±2,  and  ±4,  where  m 
indieates  azimuthal  mode.  A  detailed  deseription  of  the  azimuthal  modes  is  in  Kim  et  al.  [2009a]. 
Although  experiments  are  eondueted  for  all  these  modes,  more  extensive  results  for  m  =  0,  1,  and  ±1  will 
be  presented  sinee  these  modes  were  representative  in  Mj  =  0.9  jets  [Kim  et  al.  2009a].  The  foreing 
Strouhal  number  (SIdf  =  fFE)/Ue,  fF  is  foreing  frequeney)  ranges  from  0.07  to  3.0,  eovering  the  jet  eolumn 
mode  instability  and  the  lower  end  range  of  the  initial  shear  layer  instability.  The  jet  exit  veloeity  is  used 
in  ealeulating  the  foreing  Strouhal  numbers  for  all  jet  Maeh  numbers,  and  its  value  varies  due  to  variation 
of  the  stagnation  temperature. 

2.4  Flow  Field  Measurements 

The  jet  veloeity  field  is  measured  by  a  LaVision  PIV  system  using  either  one  or  two  eameras  with 
2048x2048  pixel  resolution.  A  Speetra  Physies  Model  SP-400  dual-head  Nd:YAG  laser  is  used  for  the 
light  souree.  The  eameras  and  laser  are  synehronized  by  a  timing  unit  housed  in  a  dual-proeessor  PC.  The 
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setup  for  the  PIV  and  the  flow  visualization  is  the  same  and  is  depieted  in  Fig.  2.4.  The  spatial  resolution 
of  the  veloeity  veetors  depends  on  the  field  of  view,  and  the  number  of  pixels  used.  For  the  most  of 
streamwise  veloeity  field  measurements,  the  spatial  resolution  is  about  2.2  -  2.5  mm. 
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Fig.  2.4  Schematic  of  the  jet  and  the  optical  diagnostics  set  up  at  GDTL.  Y-coordinate  is  normal  to 

the  plane. 

The  eold  jet  plume  is  seeded  with  Di-Ethyl -Hexyl-Sebaeat  (DEHS)  liquid  droplets  atomized  by  a 
four  jet  LaVision  atomizer.  However,  it  is  not  possible  to  seed  the  heated  jet  plume  with  liquid  droplets 
as  the  high  temperatures  of  the  jet  evaporate  the  droplets,  and  so  it  is  neeessary  to  use  an  alternative 
seeding  teehnique.  The  seeding  used  in  heated  jets  is  aluminum  oxide  partieulates  suspended  in  ethanol,  a 
teehnique  developed  and  used  by  Wernet  &  Wernet  [1994].  The  partieles  have  a  mean  diameter  of  0.6  pm 
in  a  0.02%  eoneentration  by  weight.  A  stable  suspension  is  possible  beeause  the  aluminum  oxide  and  the 
ethanol  have  signifieantly  different  pH  values.  The  ethanol,  or  some  substituted  liquid,  sometimes  has  to 
be  pH  adjusted  in  order  to  obtain  the  proper  relationship  to  the  partieles.  This  pH  imbalanee  ereates  a 
slight  eleetrie  repulsion  between  individual  partieles.  The  partieles,  whieh  on  their  own  tend  to 
agglomerate,  are  plaeed,  in  high  eoneentration,  into  ethanol  and  the  mixture  is  sonieated.  The  ultrasonie 
waves  break  apart  the  agglomerated  partieles  ereating  the  suspension.  The  eoneentrated  suspension  may 
then  be  diluted  down  to  the  desired  level  by  adding  more  ethanol.  When  done  properly,  this  suspension  is 
stable  almost  indefinitely.  In  the  experiments  eondueted  for  this  paper,  suspensions  were  sometimes  left 
to  sit  for  months  and  were  still  usable. 

A  38.1  em  (15”)  duet  is  plaeed  upstream  of  the  jet  exit  to  generate  a  eo-flow.  The  eo-flow  is 
generated  by  ehanneling  part  of  the  entrained  air  into  the  jet  through  the  duet  without  using  any  fans  or 
blowers.  The  eo-flow  is  seeded  by  a  fogger  to  avoid  spurious  veloeity  veetors  in  the  entrained  air  region. 
The  average  droplet  size  is  about  0.25  and  0.7  pm  for  the  jet  flow  and  eo-flow,  respeetively.  The 
turbulenee  statisties  were  eonverged  using  600  to  650  image  pairs  [Kim  et  al.  2009a&b].  Thus,  about  700 
image  pairs  are  used  for  all  the  statisties  reported  in  this  paper.  The  uneertainty  in  the  PIV  measurements 
is  related  to  many  parameters  sueh  as  the  partiele  size  and  density,  and  turbulenee  seales  of  interest. 
Within  5%  deviation  from  the  aetual  turbulenee  intensity,  the  seeded  partieles  traee  the  flow  up  to  20  and 
70  kHz  of  turbulenee  fluetuations  frequeney  for  0.7  and  0.25  pm  partieles  respeetively  [Melling  1997]. 
Based  on  this  ealeulation,  the  uneertainty  of  turbulenee  intensity  is  about  5%  up  to  a  Strouhal  number  of 
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1.33.  However,  the  uneertainty  level  for  the  mean  and  turbulenee  statisties  was  within  ±3%  and  ±15%, 
respeetively,  based  on  the  repeatability  measurements  for  the  baseline  jet.  In  the  shoek-eontaining 
imperfeetly  expanded  jets,  the  partieles  lag  behind  the  aetual  flow  speed  in  regions  near  the  shoeks. 
Melling  [1997]  showed  that  a  0.25  pm  partiele  passing  through  an  oblique  shoek  wave  (upstream  and 
downstream  Maeh  numbers  are  1.5  and  1.15,  respeetively  needed  about  0.5  mm  before  it  was  within  5% 
of  the  downstream  veloeity.  In  the  present  researeh,  the  shoek  is  not  quite  as  strong  as  Melling’ s  example 
so  an  estimated  distanee  required  for  the  partiele  to  reaeh  the  downstream  veloeity  is  about  0.2  mm.  Thus, 
the  uneertainty  of  the  present  PIV  measurements  is  not  affeeted  by  partiele  lag  for  most  of  the  flow  field 
measured. 
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3.  RESULTS 

The  three  major  eontrol  parameters  are  duty  eyele,  foreing  frequeney  (or  Strouhal  number),  and 
azimuthal  mode.  The  effeets  of  these  parameters  will  be  investigated  in  eold  Mj  =  0.9,  1.2,  1.3,  and  1.4 
jets.  The  Maeh  0.9  jet  was  used  to  further  develop  both  eontrol  and  measurements  tools  before  moving 
on  to  supersonie  jets.  The  measures  for  the  jet  spreading  used  in  the  researeh  are  eenterline  Maeh  number 
and  the  jet  width  at  half  maximum  (the  jet  width  based  on  the  half  veloeity  of  the  loeal  eenterline 
veloeity).  Large-seale  struetures  generated  by  foreing  are  visualized  and  their  role  in  jet  mixing 
enhaneement  are  also  diseussed.  Four  main  results  to  be  diseussed  are: 

a)  Mj  =  0.9  unheated  subsonie  jet  (extensive  foreing  frequeneies  and  azimuthal  modes.  See.  3.1) 

b)  Mj  =  1.2,  1.3,  and  1.4  unheated  supersonie  jet  (limited  azimuthal  modes.  See.  3.2) 

e)  Mj  =  0.9  heated  subsonie  jet  (limited  azimuthal  modes  at  single  foreing  frequeney.  See.  3.3) 

d)  Mj  =  1.65  unheated  subsonie  jet  (only  for  one  mode  and  several  foreing  frequeney.  See.  3.4) 


As  diseussed  in  See.  1.3,  eaeh  seetion  has  unique  objeetives  as  follows: 

3.1)  To  elarify  the  effeets  of  duty  eyele,  foreing  frequeney,  and  azimuthal  modes  in  Mj  =  0.9  unheated 
jets  and  to  find  the  optimum  duty  eyele  for  eaeh  foreing  frequeney. 

3.2)  To  determine  the  eontrol  authority  of  plasma  aetuators  in  unheated  supersonie  jets  of  Mj  =  1.2, 
1.3,  and  1.4  for  limited  azimuthal  modes  based  the  results  in  the  Mj  0.9  jet. 

3.3)  To  find  the  effeets  of  the  stagnation  temperature  on  the  effeetiveness  of  plasma  aetuators  in  Mj  = 
0.9  heated  jets  for  very  limited  azimuthal  mode  and  fixed  foreing  Strouhal  number  of  0.3. 

3.4)  To  determine  the  eontrol  authority  in  a  high  Maeh  number  jet  of  Mj  =  1.65  perfeetly-expanded 
unheated  jet  only  for  m  =  ±1  and  several  Stop’s. 


The  extensive  results  in  Mj  =  0.9  and  1.3  jets  serve  as  a  referenee  for  generating  a  test  matrix  for  other 
experiments. 


3.1  Mach  0.9  Cold  Subsonic  Jet 

Detailed  two-eomponent  partiele  image  veloeimetry  measurements  on  a  streamwise  plane  passing 
through  the  jet  eenterline  are  used  to  explore  the  effeets  of  foreing  Strouhal  number  and  azimuthal  mode 
on  the  Maeh  0.9  jet.  The  overall  performanee  of  the  plasma  aetuators  at  eaeh  Strouhal  number  and  mode 
are  diseussed  by  using  the  average  veloeity  images  and  turbulenee  statisties.  Then,  eonditionally- 
averaged  veloeity  eomponents  superimposed  on  eonditionally-averaged  streamlines  are  be  used  to 
investigate  the  dynamies  of  vortiees  or  large-seale  struetures  and  their  roles  in  the  jet  development. 


3.1.1  Effects  of  Duty  Cycle  of  Input  Signal 

The  effeet  of  duty  eyele,  pereentage  of  are  duration  to  pulsation  period,  was  investigated  at 
various  Stores  for  m  =  ±1.  However,  only  results  for  Stop  =  0.32,  near  the  jet  eolumn  instability,  are 
shown  in  Fig.  3.1.1.  The  figure  shows  the  Maeh  number  deeay  along  the  jet  eenterline  up  to  x/D  =  1 1  at 
duty  eyeles  ranging  from  3%  to  50%.  The  eenterline  Maeh  number  is  an  indireet  indieator  of  jet  spreading 
and  the  deeay  of  Maeh  number  is  faster  for  inereased  jet  spreading  in  general.  As  the  duty  eyele  deereases. 
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the  centerline  Mach  number  decays  faster,  implying  faster  jet  spreading.  When  the  duty  cycle  is  further 
decreased,  there  is  occasional  misfire  in  plasma  actuators.  Thus,  the  optimal  duty  cycle  at  a  forcing 
frequency  is  determined  at  smallest  duty  cycle  which  does  not  cause  misfire  in  actuators.  The  following 
equation  shows  the  relation  between  duty  cycle  (t  is  not  in  percentage  but  in  fraction,  so  t  =  0  -  1)  and 
forcing  frequency  (fp). 


J [0.6(/^  / 1 000)+2] / 1 00  for  <30 kHz 

I [0.26(/^  / 1 000)+ 1 1 .42]  / 1 00  for  f^  >  30kHz. 


In  all  following  experiments,  the  duty  cycle  is  determined  by  Eq.  1.  Then  the  pulse  train  g(t)  for  each 
actuator  is  given  as: 


0  if  0<t< 


1  if 

0  if 


Ijfp 


Ijfp 

m^a 


2+f  fp 


+  —  <?  <1 
2+f  fp 


(2) 


where  (l)^=Ti{i-l)/4  is  the  channel  phase  corresponding  to  the  azimuthal  location  of  the  actuator  in  the 

nozzle  extension  (see  Fig.  2.1  or  2.3),  i=l,  2,  N  is  the  channel  number,  N  is  the  number  of  channels 
(7V=8  in  the  present  work),  and  m  is  the  phase  shift  parameter  (azimuthal  mode  number).  For  m  =  0,  all 
plasma  actuators  operate  in  phase,  producing  axisymmetric  forcing  mode;  m=\  produces  sequential 
operation  in  the  first  helical  mode;  and  m=±\  generates  two  counter-rotating  helical  modes,  superposition 
of  which  produces  a  “flapping”  mode.  A  more  detailed  description  of  the  input  pulse  train  is  given  in 
Samimy  et  al.  [2007a]  and  Kim  et  al.  [2009a]. 


Fig.  3.1.1  Centerline  Mach  number  at  Stop  =  0.32  (3.5  kHz)  and  at  m  =  ±1. 
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3.1.2  Effects  of  Forcing  Strouhal  Number  on  Overall  Jet  Spreading 

The  jet  was  foreed  at  foreing  Strouhal  numbers  (Stop)  ranging  from  0.09  to  3.08  for  all  available 
azimuthal  modes.  However,  only  results  for  m  =  0,  1,3  and  ±1  are  shown  in  this  seetion.  The  eenterline 
veloeity  deeay  has  been  widely  used  as  a  measure  of  the  jet  spreading  or  mixing  with  the  ambient  air,  as 
often  it  is  the  only  available  results.  The  eenterline  veloeity  deeay  inereases  as  the  jet  spread  inereases  in 
most  eases.  In  the  present  researeh,  growth  of  the  full  width  at  half  maximum  (FWHM)  of  the  jet 
eenterline  veloeity  is  also  used  as  a  measure  for  the  jet  spreading  or  mixing. 


Fig.  3.1.2  Jet  width  development  with  downstream  locations  (a),  centerline  Mach  number  decay 
(b),  and  centerline  TKE  (c)  at  m  =  0. 

For  m  =  0  mode,  the  jet  width  (ealeulated  by  using  FWHM  and  referred  as  the  jet  width 
heneeforth),  is  shown  in  Fig.  3.1.2  along  with  the  jet  eenterline  Maeh  number  and  the  measured  2- 
eomponent  turbulenee  kinetie  energy  (from  here  on  TKE)  profiles.  The  trend  of  jet  width  with  the  Stop 
does  not  mateh  the  eenterline  Maeh  number  deeay  trend  at  this  mode.  However,  the  trend  of  the  eenterline 
Maeh  number  deeay  is  quite  similar  to  that  of  the  eenterline  TKE  development.  As  will  be  shown  later, 
the  mismateh  of  these  trends  is  observed  only  at  the  axisymmetrie  mode  (m  =  0).  At  this  point,  it  is  not 
elear  why  they  are  so  different  at  this  mode. 

The  jet  width  shows  the  effeets  of  foreing  on  the  jet  spreading  for  the  entire  streamwise 
measurement  extent,  whereas  the  eenterline  Maeh  number  deeay  provides  the  overall  mixing/spreading 
effeet  only  beyond  the  potential  eore.  When  the  jet  is  foreed  at  a  Strouhal  number  near  or  below  0.36,  the 
jet  width  grows  almost  linearly  in  downstream  direetion  for  the  entire  measured  x/D.  The  best 
spreading/mixing  is  observed  at  Stop  =0.27,  where  the  jet  width  is  inereased  by  about  20%  at  x/D  =  9.  The 
jet  width  grew  faster  and  saturated  eloser  to  the  nozzle  exit  when  the  jet  was  foreed  at  higher  Strouhal 
numbers  resulting  in  signifieant  mixing  reduetion  by  x/D  =  9. 

As  the  foreing  Strouhal  number  is  inereased,  the  saturation  oeeurred  earlier  sinee  the  generated 
struetures  are  smaller  and  thus  their  life  span  is  shorter.  At  a  Stop  equal  or  greater  than  1.0,  the  jet  width  is 
redueed  upstream  of  x/D  =  2  and  remained  unehanged  up  to  x/D  =  6  (whieh  is  approximately  the  end  of 
the  potential  eore).  The  jet  spreading  is  redueed  by  foreing  at  a  higher  Sto  most  likely  by  suppressing  the 
formation  and/or  development  of  large-seale  struetures,  whieh  play  an  important  role  in  entrainment  and 
jet  spreading/mixing. 

Figure  3.1.3  shows  the  development  of  jet  width  and  the  eenterline  Maeh  number  at  m  =  1  mode. 
For  this  mode,  unlike  for  the  m  =  0  mode,  the  trend  of  the  eenterline  Maeh  number  is  similar  to  that  of  the 
jet  width  development  downstream  of  x/D  of  about  5.  Although  not  shown,  the  trend  of  the  eenterline 
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TKE  matches  with  these  two  trends  as  well.  As  mentioned  earlier,  the  three  trends  are  very  similar  for  all 
the  azimuthal  modes  tested  except  for  m  =  0. 


(a)  Jet  width  (b)  Centerline  Mach  number 

Fig.  3.1.3  Growth  of  the  jet  width  with  downstream  locations  (a)  and  centerline  Mach  number 

decay  (b)  at  m  =  1. 

For  m=l,  the  growth  of  the  jet  width  is  maximum  when  forcing  at  Stop  =  0.18,  and  only  a  slightly 
lower  growth  is  observed  at  Stop  =  0.27  and  0.32.  For  these  three  forcing  cases,  as  well  as  others  with 
relatively  strong  jet  spreading,  the  initial  linear  growth  of  the  jet  width  stops  near  the  end  of  potential  core 
before  starting  another  linear  growth  region,  but  with  much  larger  slope.  At  a  Stop  of  0.36  (not  shown 
here)  and  higher,  there  was  a  region  of  almost  zero  growth  rate  as  clearly  observed  at  a  Stop  of  1.09.  The 
significantly  reduced  growth  near  the  end  of  the  jet  potential  core  is  due  to  the  cross-centerline  interaction 
of  large  scale  structures  in  this  region.  This  will  be  further  discussed  later. 

At  Stop’s  greater  than  1.0,  the  initial  growth  rate  of  the  jet  width  is  increased  over  that  of  the 
baseline  (Fig.  3.1.3a),  but  decreased  for  Stop’s  less  than  0.7.  The  saturation  in  the  jet  width  growth  occurs 
earlier  and  the  saturated  jet  width  also  decreases  as  the  Stop  increases,  which  is  similar  to  what  was  found 
at  m  =  0  mode.  In  the  near  field,  the  growth  and  saturation  of  the  jet  width  with  the  Stop  is  consistent  with 
the  pressure  perturbation  level  measurements  along  the  nozzle  lip  line  using  a  single  actuator  [Samimy  et 
al.  2007].  The  jet  width  remains  about  the  same  after  it  saturated,  and  the  jet  growth  is  suppressed  at  high 
Stop’s  as  was  true  for  the  m  =  0  mode.  As  will  be  discussed  in  a  later  section,  the  quick  jet  width 
saturation  with  increasing  Stop  is  due  to  the  smaller  structures  and  their  shorter  life  span  when  forcing  the 
jet  at  higher  frequencies. 

For  the  first  combined  mode  (m  =  ±1),  often  referred  to  as  the  flapping  mode,  the  jet  width 
growth  on  the  flapping  plane  is  shown  in  Fig.  3.1.4a.  The  trend  of  the  jet  width  growth  is  consistent  with 
the  centerline  Mach  number  decay  as  shown  in  Fig.  3.1.4b.  Note  that  the  ordinate  scale  in  Fig.  3.1.4a  is 
over  twice  that  for  the  other  modes  presented  earlier.  In  the  upstream  region,  the  trend  of  the  jet  width 
growth  with  Stop  is  similar  to  that  for  the  m  =  1  mode.  Similar  to  the  results  for  the  other  modes  presented, 
the  growth  in  jet  width  is  suppressed  when  the  jet  is  forced  at  higher  Stop’s,  as  can  be  seen  at  3.08.  The  jet 
width  is  increased  significantly  near  Stop  =  0.3,  and  the  maximum  growth  is  at  Stop  ^  0.27,  approximately 
3  times  that  of  the  baseline  at  x/D  =  9.  At  low  Stop’s,  the  jet  width  grows  monotonically  up  to  x/D  =  4, 
which  is  near  the  end  of  the  potential  core.  When  Stop  is  increased  or  decreased  from  0.27,  the  jet 
spreading  is  decreased  very  rapidly.  At  Stop  of  0.09  and  0.73,  the  jet  width  increases  monotonically 
(almost  linearly)  with  downstream  location.  The  growth  in  jet  width  near  Stop  ^  0.3  starts  to  increase 
exponentially  at  x/D  =  4,  where  the  potential  core  ends  as  seen  in  Fig.  3.1.4b.  As  with  the  other  azimuthal 
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modes,  the  enhaneed  growth  resulted  in  a  reduetion  in  the  jet  potential  eore.  The  potential  eore  length  was 
redueed  to  x/D  =  4  for  the  best  mixing  ease  from  x/D  =  6  for  the  baseline  ease  (Fig.  3.1 .4b). 
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01^3456739 

i/D 

(a)  Jet  width 


Fig.  3.1.4  Development  of  jet  width  and  centerline  Mach  number  for  m  =  ±1  mode  on  the  flapping 
plane.  Note  that  the  ordinate  scale  in  (a)  is  over  twice  the  others  in  Figs.  3.1.2  &  3.1.3. 


As  was  presented  and  will  be  further  diseussed  below,  the  most  effeetive  foreing  was  at  m  =  ±1 
mode  as  far  as  the  jet  spreading  is  eoneerned.  For  this  mode,  the  average  streamwise  veloeity  eontours,  on 
the  flapping  plane  of  the  jet,  are  shown  in  Fig.  3.1.5  at  Stop  =  0.18,  0.27,  0.72,  1.08,  and  3.08.  The  jet  exit 
veloeity  is  about  280  m/s  and  varies  slightly  depending  on  the  jet  stagnation  temperature.  In  the  figure, 
the  low-speed  baekground  is  for  the  eo-fiow  where  the  veloeity  is  less  than  3  m/s  (about  1%  of  the  jet  exit 
veloeity)  and  is  not  expeeted  to  affeet  the  jet  development  signifieantly. 


(d)  Stop  =  0.72  (e)  Stop  =  1.08  (f)  Stop  =  3.08 


0  50  100  150  200 


Fig.  3.1.5  Average  streamwise  velocity  fleld  for  various  forcing  Strouhal  numbers  at  m  =  ±1.  The 
velocity  scale  is  in  m/s  and  the  exit  jet  velocity  is  about  280  m/s. 
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As  Stop  increases,  the  jet  spreading  angle  increases  until  it  reaches  a  maximum  around  Stop  =  0.3, 
beyond  which  it  decreases  slowly.  The  jet  spreading  is  almost  the  same  as  that  of  the  baseline  jet  when 
forcing  at  Stop’s  greater  than  about  1.0.  The  results  suggest  that  forcing  around  a  Sto  of  0.3  is  the  most 
effective  in  spreading  the  jet,  which  is  the  jet  preferred  Sto  at  this  mode.  At  this  Stop  (0.3),  the  jet 
potential  core  length  is  significantly  shortened  as  shown  in  Fig.  3.1.5(c).  This  can  be  more  readily  seen  in 
Fig.  3.1.4b,  which  shows  the  centerline  Mach  number  distributions  with  various  forcing  Strouhal  numbers. 
The  centerline  Mach  number  decays  the  fastest  at  Stop  =  0.27.  This  is  consistent  with  the  highest  jet 
spreading,  shown  in  Fig.  3.1.5(c).  As  the  Stop  increases  beyond  about  0.3,  the  potential  core  increases, 
and  becomes  almost  the  same  as  that  of  the  baseline  jet  for  Stop’s  greater  than  about  1.0. 

The  jet  response  to  the  forcing  with  plasma  actuators  is  dependent  on  the  forcing  frequency,  duty 
cycle,  and  azimuthal  mode.  The  effects  of  duty  cycle  were  shown  in  Sec.  3.1.1,  and  the  duty  cycle  was 
predetermined  by  Eq.  1  for  each  forcing  frequency.  The  optimum  forcing  Strouhal  number,  where  the  jet 
spreading  is  maximum,  depends  on  the  forcing  azimuthal  mode  as  can  be  inferred  from  the  results 
presented  so  far.  Table  2  shows  the  optimal  forcing  Strouhal  number  for  each  mode  based  on  the  jet  width 
downstream  of  the  jet  potential  core.  The  optimum  Stop  is  about  0.3  except  for  azimuthal  modes  1  and  3. 
However,  the  Strouhal  number  for  the  second  best  is  0.27  for  m  =  1  mode  and  the  jet  width  growth  at  this 
value  is  very  close  to  the  optimum  value.  The  only  exception  is  m  =  3  mode,  which  has  a  maximum  jet 
width  growth  at  Stop  =  0.09  (1/3  of  other  cases).  The  jet  preferred  Strouhal  number,  reported  in  the 
literature  for  the  past  10-20  years  varies  from  0.2  to  0.6.  Thus,  the  optimum  Stop  for  each  mode  is  within 
the  range  in  the  literature  except  for  m  =  3. 


Table  2  Optimum  forcing  Strouhal  number  for  each  mode  based  on  the  jet  width  downstream  of  jet 
potential  core  region.  _ _ _ _ _ _ 


Mode 

0 

1 

2 

3 

±1 

±2 

±4 

Stop 

0.27 

0.18 

0.36 

0.09 

0.27 

0.32 

0.32 

3.1.3  Effects  of  Azimuthal  Mode 

As  was  presented  earlier,  the  jet  preferred  Strouhal  number,  where  the  jet  spreading/mixing  is 
maximized,  varied  for  different  azimuthal  modes,  but  remained  close  to  0.3  for  most  modes  except  for  m 
=  3,  which  was  0.09.  The  average  streamwise  velocity  contours  are  shown  in  Fig.  3.1.6  for  m  =  0-3,  ±1, 
±2,  and  ±4  modes  at  Stop  corresponding  to  those  in  Table  2.  For  all  modes,  the  jet  spreading  is 
significantly  enhanced  compared  to  that  of  the  baseline  jet  (Fig.  3.1.6e),  with  the  largest  enhancement  at 
m  =  ±1  mode  (Fig.  3.1.6f).  A  shorter  visual  potential  core  length  is  observed  for  greater  jet  spreading. 
This  is  more  readily  observed  in  the  centerline  Mach  number  profiles  shown  in  Fig.  3.1.7a.  The  length  of 
the  potential  core  for  m  =  1  and  ±2  modes  is  reduced  by  a  similar  amount.  The  least  spreading  is  observed 
when  exciting  m  =  ±4.  As  was  discussed  earlier,  the  mean  velocity  contours  and  the  centerline  Mach 
number  only  show  the  effects  of  forcing  on  the  potential  core  length  and  the  overall  spreading 
downstream  of  the  potential  core. 

The  jet  width  and  equivalent  width  development  are  shown  in  Figs.  3.1.7b  and  3.1.7,  which  show 
the  effects  of  forcing  over  the  entire  region  of  the  jet  flow  from  the  nozzle  exit  to  x/D  =  9.  Figure  3.1.7b 
shows  the  jet  width  calculated  from  Fig.  3.1.6,  and  thus  only  shows  the  jet  width  on  the  PIV  plane  in  this 
figure.  Since  the  jet  cross-section  at  a  far  downstream  location  is  nearly  axisymmetric.  Fig.  3.1.7b  is 
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appropriate  for  evaluating  the  overall  spreading  exeept  for  m  =  ±1  mode.  For  m  =  ±1  mode,  the  jet 
spreads  signifieantly  more  in  the  flapping  plane  eompared  to  the  non-flapping  plane.  The  eross-seetion  is 
assumed  to  be  approximately  elliptie  and  thus  the  equivalent  jet  width  for  m  =  ±1  is  ealeulated  from  the 
square  root  of  the  multiplieation  of  jet  widths  on  the  flapping  and  non-flapping  planes,  to  take  into 
aeeount  the  highly  non-axisymmetrie  spreading.  The  results  for  all  the  modes  are  shown  in  Fig.  3.1.7e. 
The  trend  of  equivalent  jet  width  is  eonsistent  with  eenterline  Maeh  number  exeept  for  m  =  0  mode.  This 
is  due  to  the  relatively  slow  deeay  of  the  eenterline  Maeh  number  for  m  =  0  mode  as  will  be  diseussed 
later. 


(a)  m  =  0 


(b)  m  =  1 


(e)  m  =  2 


(d)  m  =  3 
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(e)  Baseline 


(f)  m  =  ±1 


(g)  m  =  ±2 
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(h)  m  =  ±4 


Fig.  3.1.6  Average  streamwise  velocity  contours  for  various  azimuthal  modes  at  the  Stop’s  shown  in 

Table  2. 
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Fig.  3.1.7  Development  of  centerline  Mach  number  and  jet  width  and  equivalent  width  at  various 
azimuthal  modes  at  Stores  shown  in  Table  2. 


The  equivalent  jet  width  also  shows  that  the  maximum  and  minimum  jet  spreading  oeeurs  at  m  = 
±1  and  m  =  ±4  modes,  respeetively.  For  m  =  2,  ±1,  and  ±4  modes,  the  jet  width  grows  signifieantly  over 
the  other  modes  in  the  initial  region  or  near  field.  However,  the  growth  rate  of  the  jet  width  for  m  =  2  and 
±4  deereases  near  the  end  of  the  potential  eore  and  thus  this  lead  to  a  limited  inerease  in  jet  widths  by  x/D 
=  9  (Fig.  3.1.7e).  For  m  =  ±2  mode,  the  jet  width  remains  saturated  from  x/D  =  2.5  until  x/D  =  5 
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(approximate  end  of  the  potential  eore),  and  then  experienees  rapid  growth  further  downstream.  At  higher 
azimuthal  modes  of  m  =  3  and  ±4,  the  jet  eenterline  Maeh  number  deeay  is  slower  and  the  inerease  in  the 
jet  width  at  a  far  downstream  loeation  is  less  than  those  of  other  modes.  For  m  =  ±2,  the  eenterline  Maeh 
number  deeay  is  eomparable  to  that  of  m  =  ±1  although  the  streamwise  veloeity  eontour  is  not 
eomparable.  In  addition  to  the  jet  spreading,  it  seems  that  the  vortex- vortex  interaetion  aeross  the  jet 
eolumn  is  another  faetor  related  to  the  eenterline  Maeh  number  deeay.  This  will  be  further  diseussed  later 
by  using  large-seale  struetures  and  their  dynamie  interaetion. 


3.1.4  Effects  Forcing  on  Turbulence  Statistics 

As  presented  in  the  previous  seetion,  ehanges  in  the  mean  flow  eharaeteristies  depend 
signifieantly  on  the  Stop  and  foreing  azimuthal  mode.  In  this  seetion,  the  effeets  of  foreing  on  the 
turbulenee  statisties  are  explored  along  the  jet  eenterline.  The  normalized  two-dimensional  turbulent 
kinetie  energy  and  anisotropy  ratio  (Qy/au,  and  Gy  are  RMS  of  x-  and  y-eomponent  veloeity 
fluetuations,  respeetively)  along  the  jet  eenterline  are  shown  in  Fig.  3.1.8  for  m  =  0,  1,  and  3  modes. 


(d)  Anisotropy  for  m  =  0 


(e)  Anisotropy  for  m  =  1 


(f)  Anisotropy  for  m  =  3 


Fig.  3.1.8  Normalized  TKE  and  anisotropy  for  m  =  0, 1,  and  3  modes  along  the  jet  centerline. 


For  m  =  0  mode,  the  TKE  is  inereased  for  all  Stop’s  with  the  maximum  amplifieation  oeeurring  at 
Stop  =  0.18,  whieh  eorresponds  to  the  ease  where  the  eenterline  Maeh  number  deeayed  the  fastest  (Fig. 
3.1.2b).  The  TKE  is  more  amplified  at  Stop  =  0.32  and  0.36  (not  shown  here)  up  to  around  the  end  of  the 
potential  eore  (x/D~  5).  Up  to  x/D  =  1  to  1.5,  the  anisotropie  ratio  is  1,  as  expeeted,  as  the  fluetuations 
are  due  to  random  noise  and  measurement  errors.  However,  further  downstream,  the  ratio  is  substantially 
less  than  one  at  StDP<0.73,  and  substantially  larger  than  one  for  the  baseline  jet  and  at  StDP>1.09  for  m=0 
ease  (Fig.  3.1.8d).  At  lower  Stop’s,  the  streamwise  (or  x-eomponent)  turbulenee  intensity  was  amplified 
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more  than  the  eross-stream  (y-eomponent)  turbulenee  intensity.  At  higher  Stop’s,  the  seales  of  the 
generated  struetures  are  so  small  that  their  influenee  on  the  eenterline  is  expeeted  to  be  very  limited,  and 
thus  the  foreing  effeets  are  expeeted  to  be  negligible.  On  the  other  hand,  the  seales  of  the  indueed 
struetures  are  mueh  larger  at  lower  Stop’s,  inereasing  interaetions  aeross  the  jet  eenterline  and  thus 
signifieantly  altering  the  anisotropy  ratio.  The  anisotropy  ratio  was  redueed  for  all  eases  downstream  of 
the  potential  eore. 

For  m  =  1,  a  signifieant  amplifieation  in  TKE  was  observed  over  a  wide  range  of  Stop’s  from  0.18 
to  0.36  (not  shown  here).  The  TKE  is  saturated  around  x/D  =  7  for  these  foreing  Strouhal  numbers. 
Contrary  to  the  m  =  0  ease,  the  anisotropy  in  Fig.  3.1.8e  implies  that  the  eross-stream  veloeity  fluetuations 
are  dominant  over  the  streamwise  veloeity  fluetuations  for  Stop’s  =  0.18-0.36  (0.36  ease  is  not  shown)  in 
the  potential  eore  region.  As  will  be  diseussed  later,  the  differenee  is  assoeiated  with  the  symmetrie  or 
asymmetrie  nature  of  the  large-seale  struetures  aeross  the  jet  diameter.  The  turbulenee  field  beeomes 
more  isotropie  downstream  of  the  potential  eore.  Figure  3.1.8e  shows  that  the  anisotropy  is  saturated  near 
the  end  of  potential  eore.  As  the  Stop  is  inereased,  the  amplifieation  level  deereases  and  the  anisotropy  is 
about  the  same  level  as  the  baseline  jet  ease  for  Stop‘s  equal  or  greater  than  1.8  (not  shown  here). 

For  m  =  3,  the  amplifieation  in  TKE  is  moderate  exeept  for  Stop  =  0.09  and  0.18,  where  the  jet 
spreading  was  maximum.  At  a  low  Stop,  the  anisotropy  inereased  near  the  end  of  the  potential  eore 
similar  to  that  for  m  =  1,  but  the  inerease  is  moderate.  Interestingly,  the  TKE  level  was  redueed  when  the 
jet  was  foreed  at  a  Stop  greater  than  1.0.  The  reduetion  in  TKE  at  high  Stop’s  is  related  to  the  broadband 
noise  suppression  seen  for  this  mode  [Samimy  et  al.  2007a,  Bridges  and  Wernet  2002]. 

The  development  of  TKE  and  anisotropy  along  the  jet  eenterline  for  m  =  ±1  are  shown  in  Fig. 
3.1.9.  A  dramatie  inerease  in  TKE  was  observed  at  low  Stop’s.  The  anisotropy  also  inereased  signifieantly 
at  low  Stop’s,  but  the  trend  showed  some  disparity  from  that  of  TKE.  The  inerease  in  anisotropy  means 
that  the  veloeity  fluetuations  in  the  eross-stream  direetion  are  dominant  over  those  in  the  streamwise 
direetion,  as  was  also  seen  for  m  =  1 .  Sinee  the  flapping  plane  is  in  the  y-direetion,  the  domination  of  the 
eross-stream  veloeity  fluetuations  is  expeeted.  As  the  Stop  is  inereased,  the  development  of  TKE  and 
isotropy  approaeh  the  levels  of  the  baseline  jet. 


(a)  TKE  (e)  Anisotropy 

Fig.  3.1.9  Normalized  TKE  and  anisotropy  for  m  =  ±1  mode  along  the  jet  centerline. 

Figure  3.1.10  shows  the  development  of  TKE  and  anisotropy  at  the  modes  and  Stop’s,  where  the 
jet  width  was  maximized  as  shown  in  Table  2.  The  interaetion  between  the  generated  struetures  aeross  the 
jet  eenterline  is  expeeted  to  be  notieeable  when  the  seale  of  struetures  is  large  and  eomparable  to  the  jet 
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diameter.  The  onset  of  where  the  TKE  starts  to  inerease  over  the  baekground  level  appears  to  be 
approximately  the  same  for  all  azimuthal  modes,  but  the  amplifieation  rates  are  signifieantly  different  for 
eaeh  mode.  For  m  =  0,  the  TKE  developed  monotonieally  with  downstream  loeation,  and  did  not  show 
any  saturation  in  the  measured  range.  Vortex  rings  were  generated  by  the  plasma  aetuators  as  will  be 
shown  later,  and  the  seale/size  of  the  rings  grew  monotonieally  with  downstream  loeation,  and  thus  the 
TKE  was  also  expeeted  to  grow  in  sueh  a  fashion.  For  m  =  ±1,  the  TKE  started  to  grow  signifieantly  as 
early  as  x/D  =  2  and  saturated  around  x/D=7.  The  least  amplifieation  in  TKE  is  observed  for  m  =  ±4  mode, 
where  its  growth  rate  is  about  the  same  as  the  baseline  jet.  For  most  eases,  the  growth  in  TKE  is  saturated 
around  x/D  =  7-8. 


(a)  TKE  (b)  Anisotropy 

Fig.  3.1.10  Development  of  TKE  and  anisotropy  at  the  Stor’s  in  Table  2. 


The  development  of  anisotropy  in  Fig.  3.1.10b  is  very  interesting.  For  even  numbered  modes, 
anisotropy  is  deereased,  implying  the  streamwise  veloeity  fiuetuations  are  dominant  in  TKE.  The  ease  for 
m  =  2  seemed  different,  but  the  development  at  a  lower  SIdf  is  similar  to  other  even  numbered  modes. 
When  the  jet  was  foreed  at  odd  numbered  modes,  the  anisotropy  is  inereased  and  saturated  around  x/D  = 
3.5.  The  anisotropy  reaehed  a  minimum  for  the  even  numbered  modes  around  x/D  of  2-2.5.  These 
differenees  ean  be  explained  through  vortex  dynamies  of  the  generated  struetures,  as  will  be  presented 
and  diseussed  in  the  following  seetion. 


3.1.5  Vortex  Dynamics  and  Its  Role  in  the  Jet  Development 

The  overall  effeets  of  SIdf  and  foreing  azimuthal  mode  were  investigated  and  diseussed  in  the 
earlier  seetions  by  using  the  average  veloeity  eontours,  the  jet  width,  TKE,  and  anisotropy.  The  average 
veloeity  and  turbulenee  statisties  are  useful  in  evaluating  the  overall  effeets  of  StoF  and  foreing  azimuthal 
mode.  However,  they  do  not  reveal  details  of  flow  struetures  and  their  role  in  the  jet  development.  Thus, 
large-seale  struetures  are  extraeted  from  PIV  data  and  their  dynamies  are  diseussed  in  this  seetion. 

The  Galilean  deeomposition  is  applied  to  the  measured  veloeity  fields  to  extraet  large-seale 
struetures.  The  eonveetion  veloeity  of  large-seale  struetures  must  be  known  to  obtain  the  Galilean- 
deeomposed  veloeity  field.  Onee  the  eonveetion  veloeity  of  large-seale  struetures  is  known,  the  Galilean- 
deeomposed  veloeity  field  is  obtained  by  subtraeting  the  eonveetion  veloeity  from  the  measured  veloeity 
field  [Konstantinidis  et  al.  2005].  Thus  in  the  Galilean  deeomposition,  the  referenee  frame  moves  at  the 
eonveetion  veloeity.  A  large-seale  strueture  does  exist  and  is  visualized  in  the  Galilean-deeomposed 
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velocity  field  if  the  streamlines  make  a  closed  loop  or  have  a  spiral  pattern  [Kline  and  Robinson  1990, 
Robinson  et  al.  1989]. 
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(a)  St  =  0.32 


(b)  St  =  0.32 
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(e)  St  =1.08  (f)  St  =1.08 

Fig.  3.1.11  Conditionally-averaged  Galilean  decomposed  velocity 
fields  and  streamlines  of  streamwise  (left  column)  and  cross- 
streamwise  (right  column)  velocities  at  m  =  ±1. 


In  the  present  research,  a  conditionally-averaged  velocity  field  was  obtained  and  used  for  large- 
scale  structures  visualization.  To  get  a  conditionally-averaged  field  from  700  instantaneous  velocity  fields, 
two-dimensional  cross-correlation  was  used  to  extract  large-scale  structures.  A  correlation  window  was 
taken  in  a  given  instantaneous  velocity  field.  The  velocity  pattern  within  the  window  was  correlated  to  the 
available  700  instantaneous  velocity  fields.  All  instantaneous  images  with  a  correlation  level  greater  than 
a  threshold  value  were  ensemble  averaged  to  form  a  conditionally-averaged  image  of  the  2-D  velocity 
field.  One  could  also  use  proper  orthogonal  decomposition  for  this  purpose  as  in  Kastner  et  al.  [2008]. 
The  next  step  is  to  calculate  convection  velocity  of  large-scale  structures.  Two-dimensional  spatial- 
correlations  were  used  to  calculate  the  spacing  or  wavelength  of  large-scale  structures  for  the  forced  cases. 
The  convection  velocity  is  obtained  by  multiplying  the  large-scale  structure  spacing  by  the  forcing 
frequency  as  was  done  by  Troutt  and  McLaughlin  [1982].  The  last  step  to  obtain  a  conditionally-averaged 
Galilean-decomposed  velocity  field  is  to  subtract  the  convection  velocity  from  the  conditionally-averaged 
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velocity  field.  Then,  the  large-scale  structures  are  visualized  when  the  Galilean  streamlines  are  added,  as 
was  detailed  in  Kline  and  Robinson  [1990]  and  in  Robinson  et  al.  [1989],  and  is  shown  in  Fig.  3.1.11  for 
m=±l  at  Stop’s  of  0.32,  0.72,  and  1.08. 

In  the  figure,  the  conditionally-averaged  Galilean  streamlines  and  streamwise  velocity  magnitude 
(Fig.  3.1.11a)  and  cross-stream  velocity  magnitude  (Fig.  3.1.11b)  are  superimposed.  For  the  streamwise 
velocity,  the  dark  and  bright  regions  represent  the  slower  and  faster  velocities,  respectively.  The 
maximum  and  minimum  velocities  are  about  280  (in  the  potential  core)  and  3  m/s  (in  the  ambient  region). 
For  the  cross-stream  velocity,  the  medium  tone  background  represents  near  zero  velocity,  and  the  brighter 
and  darker  tones  indicate  positive  and  negative  velocities,  respectively.  In  Figs.  3.1.11-13,  the  background 
velocity  contours  are  added  for  the  ease  of  visualization.  The  figure  shows  the  size  and  spacing  of  vortices, 
and  the  interaction  between  not  only  the  vortices  but  also  the  vortices  and  the  jet  column.  The  large-scale 
structures  generated  by  the  forcing  are  robust  and  seem  to  be  two-dimensional  on  the  visualization  plane, 
rotating  either  clockwise  or  counter-clockwise,  for  those  in  the  bottom  shear  layer  and  top  shear  layer, 
respectively. 

At  lower  Stop’s,  the  generated  large-scale  structures  are  very  well  organized  and  their  scale  by  the 
end  of  the  potential  core  is  comparable  to  the  nozzle  exit  diameter.  The  vortices  in  the  top  and  bottom 
shear  layer  are  out  of  phase  since  the  jet  was  forced  with  the  flapping  mode  (m  =  ±1).  Downward 
velocity  is  induced  in  the  downstream  side  of  a  vortex  in  the  bottom  shear  layer  and  the  upstream  side  of  a 
vortex  in  the  top  shear  layer,  while  upward  velocity  is  induced  in  the  upstream  side  of  a  vortex  in  the 
bottom  shear  layer  and  the  downstream  side  of  a  vortex  in  the  top  shear  layer,  as  indicated  by  arrows  in 
Fig.  3.1.11b.  The  cross-stream  velocity,  induced  by  the  vortices,  appears  strong  enough  to  cause 
undulations  in  the  jet  column,  which  can  be  inferred  from  wavy  streamlines  along  the  jet  centerline.  When 
the  vortices  in  the  top  and  bottom  shear  layers  are  out-of-phase,  as  expected  for  m=±l  case,  the  upward 
and  downward  induced  velocities  are  in  the  same  cross-stream  direction  at  the  same  streamwise  location, 
as  shown  in  Fig.  3.1.1  lb.  As  a  result,  the  jet  appears  to  be  flapping  by  the  induced  velocity.  In  addition  to 
the  undulating  motion  of  the  jet  column,  the  entrainment  of  the  ambient  air  and  the  ejection  of  jet  fluid 
into  the  ambient  by  the  induced  cross-stream  velocity  increase  the  lateral  spread  of  the  jet. 

When  the  Stop  is  increased  to  0.72,  the  spacing  of  the  adjacent  vortices  and  the  scale  of  the 
generated  structures  are  significantly  reduced  (about  halved)  as  shown  in  Fig.  3.1.11c.  At  this  Strouhal 
number,  the  interaction  between  the  vortices  in  the  top  and  bottom  shear  layers  is  weaker  due  to  reduced 
scale  of  the  generated  vortices.  The  jet  column  did  not  significantly  undulate  because  of  the  smaller 
vortices  and  weaker  interactions  between  the  vortices  and  the  jet  column  at  this  Strouhal  number.  As 
shown  in  Figs.  3.1.4  and  3.1.5,  the  jet  potential  core  was  not  significantly  changed  by  the  generated 
structures  due  to  the  limited  interaction  between  top  and  bottom  shear  layers.  As  shown  in  Fig.  3.1.1  Id, 
the  induced  cross-stream  velocity  is  confined  to  the  thin  top  and  bottom  shear  layers.  Since  the  induced 
cross-stream  velocity  is  reduced,  it  is  expected  that  the  entrainment  and  ejection  of  fluid  would  also  be 
limited  at  this  Strouhal  number.  As  a  result,  the  jet  spreading  was  less  enhanced  than  at  Stop  =  0.32. 

At  higher  Stop  of  1.08,  the  generated  vortices  are  barely  identifiable  and  are  not  organized, 
resembling  those  in  the  unforced  jet.  This  is  the  reason  for  the  mean  flow  (Figs.  3.1.4  and  3.1.5)  and 
turbulence  statistics  (Fig.  3.1.9)  for  this  forced  case  to  be  similar  to  those  for  the  baseline  jet.  For  all  other 
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modes,  the  effeets  of  Stop  on  the  spaeing  and  size  of  the  generated  stmetures  are  very  similar,  and 
therefore,  the  results  are  not  presented  here. 


(e)  m  =  1  (d)  m  =  1 

Fig.  3.1.12  Conditionally-averaged  Galilean  streamlines  and 
velocity  magnitude  contours  for  streamwise  (left  column)  and 
cross-streamwise  (right  column)  velocities  for  m  =  0  and  1 
modes  at  Stop  =  0.32. 


Figure  3.1.12  shows  eonditionally-averaged  streamwise  and  eross-stream  veloeity  eontours  with 
superimposed  Galilean  streamlines  for  two  additional  modes  m  =  0  and  1  at  a  Stop  of  0.32.  As  diseussed 
earlier,  the  foreing  at  even-  and  odd-numbered  azimuthal  modes  showed  distinetly  different  turbulenee 
eharaeteristies  as  shown  in  Fig.  3.1.10b  -  henee  the  seleetion  of  these  two  representative  modes.  For  the 
axisymmetrie  mode  (m  =  0),  the  streamwise  dimensions  of  the  generated  vortiees  are  approximately  the 
same  as  those  for  the  flapping  mode  (m  =  ±1),  but  the  eross-stream  seales  are  smaller.  Although  the 
vortiees  at  this  mode  appear  to  be  as  strong  as  those  in  the  m  =  ±1  mode,  the  jet  spreading  is  not  as 
signifieant  as  seen  in  the  m  =  ±1  ease.  The  symmetry  of  the  vortiees  seems  to  be  responsible  for  the 
slower  jet  spreading.  Sinee  the  generated  stmetures  for  m  =  0  are  donut-shaped  vortex  rings,  their 
development  in  the  eross-stream  direetion  is  limited  by  this  symmetry.  When  one  part  of  the  vortex  ring 
attempts  to  grow  toward  the  jet  eenterline,  the  opposite  part  also  takes  the  same  aetion,  and  thus  the 
growth  of  the  ring  toward  the  jet  eenterline  is  limited  due  to  the  axisymmetrie  nature  of  the  ring  vortex. 
For  the  flapping  mode  (m  =±1),  the  vortex  eould  grow  toward  the  jet  eenterline  easily  by  pushing  the  jet 
eolumn  to  the  other  side  as  shown  in  Fig.  3.1.11a.  This  limited  growth  of  the  vortex  is  partially 
responsible  for  the  slower  jet  spreading  as  seen  in  Figs.  3.1.6  and  3.1.7. 

In  addition  to  the  slower  spreading  for  the  m=0  ease,  the  jet  eenterline  Maeh  number  also  deeays 
relatively  slowly  as  seen  in  Fig.  3.1.7.  It  is  eonjeetured  that  the  jet  eenterline  Maeh  number  deeay  is 
elosely  related  to  the  interaetion  between  vortiees  and  the  jet  eolumn,  whieh  affeets  the  entrainment  of 
slow-moving  fluid  into  the  jet  plume.  The  vortiees  at  this  mode  do  not  eross  the  jet  eenterline  beeause  of 
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their  symmetry  as  diseussed  earlier.  As  a  result,  the  interaetion  is  not  as  destruetive  as  in  the  asymmetrie 
eases  (odd  numbered  modes).  The  entrained  fluid,  from  the  mueh  slower  moving  surrounding  air  near  the 
trailing  region  of  a  ring  as  indieated  by  arrows  in  Fig.  3.1.12,  goes  through  aeeeleration  at  the  eenter  of 
the  ring  due  to  self-induetion.  The  self  induetion  seems  to  make  the  flow  near  the  eenterline  aeeelerated, 
whieh  eounters  the  slowing  aetion  of  the  entrained  fluid.  Thus,  the  aeeeleration  due  to  self-induetion  at 
the  jet  eenterline  and  the  symmetry  of  the  vortex  ring  may  be  partially  responsible  for  the  slower 
eenterline  deeay  when  eompared  to  the  other  modes,  as  shown  in  Fig.  3.1.7. 

Another  interesting  flnding  at  m  =0  is  that  the  jet  eenterline  veloeity  seems  to  undulate 
periodieally  with  the  downstream  loeation.  The  jet  fluid  at  the  eenter  of  a  vortex  ring  is  aeeelerated  by 
self-induetion  of  the  vortex  ring.  In  the  trailing  region  of  a  vortex  ring,  the  veloeity  is  expeeted  to 
deerease  due  to  the  entrainment  of  the  slower  moving  ambient  air  as  indieated  by  the  arrows  in  Fig. 
3.1.12a.  The  entrainment  by  the  indueed  veloeity  ean  be  dedueed  from  the  eross-stream  veloeity  eontours, 
whieh  show  vertieally  indueed  veloeity  with  opposite  sign  as  shown  in  Fig.  3.1.12b.  Thus,  two  opposite 
aetions  of  deeeleration  and  aeeeleration  take  plaee  around  a  vortex  ring.  These  eombined  effeets  are  most 
likely  responsible  for  the  periodie  ehanging  of  the  eenterline  veloeity. 

As  diseussed  earlier,  the  development  of  TKE  along  the  jet  eenterline  for  the  m  =  0  ease  was 
monotonie  (Fig.  3.1.10a).  The  amplifieation  of  the  eenterline  TKE  is  governed  by  the  interaetion  between 
the  vortiees  and  the  jet  eolumn.  For  odd-numbered  modes,  the  flapping  aetion  of  the  jet  plume  aeross  the 
jet  eenterline  would  inerease  the  turbulenee  level.  For  m  =  0  ease,  the  vortex  rings  grew  monotonieally 
and  their  identities  were  preserved  for  long  downstream  distanees.  These  seems  to  be  the  eause  for  the 
monotonie  development  of  TKE  for  m  =  0. 

For  m  =  1  ease,  the  vortieal  struetures  are  smaller  and  weaker  than  those  for  m  =  0  or  m  =  ±1 .  The 
vortex  generated  at  this  mode  is  helieal  and  thus  the  eoherenee  level  in  the  eross-stream  direetion  is 
smaller  than  m  =  0orm  =  ±l,  where  the  generated  vortiees  are  vortieal  rings  or  spanwise  struetures, 
respeetively.  The  jet  eolumn  undulation  and  the  indueed  veloeity  around  a  vortex  are  similar  to  those  for 
the  m  =  ±1  mode,  as  ean  be  seen  in  Figs.  3.1.11  -  12.  Thus,  it  is  expeeted  that  the  vortex  dynamies  at  this 
mode  to  be  similar  to  those  at  m  =  ±1  sinee  the  vortex  patterns  are  very  similar.  However,  the  interaetion 
between  vortiees  aeross  the  jet  eolumn  and  eross-stream  veloeity  indueed  by  self-induetion  are  expeeted 
to  be  slightly  weaker  than  those  at  m  =  ±1.  This  redueed  interaetion  and  indueed  veloeity  may  be 
responsible  for  the  redueed  jet  spreading  as  shown  in  Fig.  3.1.7e. 

The  development  of  anisotropy  was  diseussed  earlier  (Fig.  3.1.10b)  and  showed  the  streamwise 
veloeity  fluetuations  to  be  largest  for  the  axisymmetrie  modes  (even  numbered  modes).  As  ean  be  inferred 
from  Fig.  3.1.12,  the  eross-stream  veloeity  fluetuations  are  suppressed  beeause  of  the  symmetrie  nature  of 
the  vortiees  aeross  the  jet  eenterline  as  shown  in  Fig.  3.1.12b.  The  indueed  veloeity  on  the  upstream  side 
of  a  pair  of  vortiees  on  the  top  and  bottom  shear  layers  is  downward  and  upward,  respeetively,  as 
indieated  by  arrows  in  Fig.  3.1.12b.  Also  the  streamwise  veloeity  fluetuation  is  more  likely  amplifled  due 
to  self-  or  mutual-induetion.  This  explains  why  the  anisotropy  is  deereased  for  the  even  numbered  modes. 
On  the  other  hand,  the  indueed  eross-stream  veloeity  fluetuations  are  in  the  same  direetion  at  a  given  x/D 
for  the  odd  numbered  modes  as  observed  in  Fig.  3. 1.1 2d.  This  leads  to  more  amplifleation  in  the  eross- 
stream  veloeity  fluetuations,  and  is  responsible  for  the  inereased  anisotropy  for  the  odd  numbered  modes. 
The  level  of  interaetion  aeross  the  jet  eenterline  ean  be  higher  earlier  for  the  even  numbered  modes  than 
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for  the  odd  numbered  ones  as  ean  be  observed  in  Figs.  3.1.12a  and  d.  This  is  assoeiated  with  the  earlier  or 
later  saturation  of  anisotropy  for  the  even  or  odd  numbered  modes,  respeetively. 


As  diseussed  in  See.  3.1.3,  the  eenterline  Maeh  number  distribution  for  m  =  ±2  is  eomparable  to 
that  of  m  =  ±1  although  the  jet  spreading  is  not.  The  vortex  dynamies  ean  offer  some  elues  for  this 
differenee.  Galilean  streamlines  superimposed  on  streamwise  and  eross-streamwise  veloeity  eomponents 
are  shown  in  Fig.  3.1.13  for  two  phases.  The  bottom  images  are  180°  out  of  phase  relative  to  the  top  ones. 
The  arrows  in  Fig.  3.1.13a  indieate  the  eenter  of  two  adjaeent  vortiees.  However,  the  arrows  at  the  same 
x/D  loeation  in  Fig.  3.1.13e  point  to  a  loeation  between  two  eonseeutive  vortiees  beeause  the  two  images 
are  out  of  phase.  Although  the  vortex  pattern  is  similar  to  that  of  the  m  =  0  mode,  the  vortiees  at  this 
mode  are  expeeted  to  be  quasi-two-dimensional  sinee  the  aetuators  on  the  vertieal  and  horizontal  planes 
are  operated  out  of  phase.  This  will  be  further  diseussed  later  using  eross-stream  veloeity  fields. 


(e)  Streamwise,  1 80°  out  of  phase  (d)  Cross-streamwise,  1 80°  out 

with  the  image  in  (a)  of  phase  with  the  image  in  (b) 

Fig.  3.1.13  Conditionally-averaged  Galilean-decomposed  velocity 
components  with  superimposed  streamlines  for  m  =  ±2. 


The  following  diseussion  is  based  on  the  assumption  that  the  generated  vortiees  are  spanwise  or 
at  least  quasi-spanwise.  This  is  eonfirmed  by  taking  images  in  the  visualized  eross-seetion  (not  shown)  of 
an  ideally  expanded  Maeh  1.3  jet  at  x/D  =  4.  The  high-speed  jet  eolumn  is  squeezed  in  the  vertieal 
direetion  between  a  pair  of  horizontally  aligned  vortiees  at  the  loeation  indieated  by  the  two  arrows  in  Fig. 
3.1.13a.  The  eorresponding  jet  eross-seetion  is  shown  sehematieally  in  Fig.  3.1.14a.  In  a  half  eyele  of  the 
foreing  period,  a  vertieally  oriented  pair  of  vortiees  pass  through  the  same  downstream  loeation  (indieated 
by  the  arrows  in  Fig.  3.1.13e)  and  eauses  the  high-speed  plume  to  be  squeezed  vertieally  as  shown 
sehematieally  in  Fig.  3.1.14b).  This  alternating  squeezing  aetion  by  vertieally  and  horizontally  aligned 
pairs  of  vortiees  is  probably  responsible  for  the  relatively  fast  deeay  of  the  jet  eenterline  Maeh  number. 
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On  the  other  hand,  it  is  expeeted  that  the  mixing/spreading  performanee  of  pairs  of  vortiees  at  this  mode 
would  not  be  large,  as  the  vortiees  are  symmetrie  aeross  the  jet  eenterline  as  explained  for  the  m  =  0  ease. 


(a)  At  an  x/D  loeation  mdieated  by  an  arrow  in  Fig.  (b)  At  an  x/D  loeation  mdieated  by  an  arrow  in  Fig. 
3.1.13a  3.1.13e 

Fig.  3.1.14  Schematic  of  the  jet  cross-section  for  m  =  ±2.  The  light  gray  tone  represents  the  jet  cross- 
section  and  the  darker  ellipses  indicate  large-scale  structures. 


Relative  x/D 

Fig.  3.1.15  Spatial  correlation  along  the  lip-line  of  the 
jet  for  several  Strouhal  numbers  at  m  =  ±1  mode. 


3.1.6  Convection  Velocity  of  Large-Scale  Structures 

The  convection  velocity  of  the  generated  structures  due  to  forcing  was  calculated  by  the  method 
discussed  in  the  earlier  section.  A  correlation  window  in  an  instantaneous  velocity  field  was  selected, 
covering  approximately  75  -  125  %  of  the  streamwise  spacing  of  two  adjacent  large-scale  structures  and 
the  entire  width  of  the  shear  layer.  As  in  Samimy  et  al.  [2007],  one  can  get  the  spatial-correlation  profiles 
along  any  streamwise  line,  e.g.  the  jet  lip-line,  from  such  two  dimensional  spatial-correlations.  The 
spatial-correlation  profiles  for  several  Strouhal  numbers  are  shown  in  Fig.  3.1.15  for  m  =  ±1.  When  there 
are  periodic  structures  in  the  shear  layer,  the  spatial  correlation  is  similar  to  an  amplitude-modulated 
sinusoidal  wave.  For  the  baseline  jet,  there  is  no  periodic  motion.  When  the  jet  was  forced,  large-scale 
periodic  structures  are  generated  as  indicated  by  the  multiple  local  peaks  of  the  amplitude -modulated 
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sinusoidal  wave.  As  discussed  earlier,  the  peak-to-peak  distance  is  the  wavelength  (spacing)  of  the 
generated  structures.  The  convection  velocity  is  obtained  by  multiplying  the  calculated  wavelength  by  the 
forcing  frequency.  The  calculated  convection  velocity  for  three  forcing  modes  of  m  =  0,  1,  and  ±1  are 
shown  in  Fig.  3.1.15. 

As  the  StDF  is  increased,  the  normalized  convection  velocity  is  increased  for  all  three  modes  as 
shown  in  Fig.  3.1.16.  At  the  lowest  forcing  Strouhal  number,  the  convection  velocity  is  close  to  the 
theoretically  predicted  value  of  0.52.  It  seems  that  the  normalized  convection  velocity  converges  to  about 
0.77  at  a  St^/s  greater  than  about  0.72.  The  present  normalized  convection  velocity,  calculated  from 
instantaneous  velocity  fields,  is  slightly  larger  than  that  obtained  from  qualitative  flow  visualization 
results  in  a  Mach  1.3  jet  for  Stop’s  of  0.3-0.65  [Samimy  et  al.  2007].  Their  results  showed  a  constant  value 
of  0.67  over  this  Strouhal  number  range.  The  convection  velocity  measured  in  a  supersonic  round  jet  at 
Mach  numbers  from  1.6  to  1.86  was  from  0.7  to  0.78  [Petitjean  et  al.  2006],  which  agrees  well  with  the 
present  results  at  higher  Strouhal  number  forcing  cases.  They  used  a  two-point  space-time  correlation 
based  technique  to  calculate  convection  velocity  in  the  cold  jet.  Bridges  [2006]  also  measured  convection 
velocity  in  cold  and  heated  circular  jets  by  using  a  dual-PIV  system.  The  normalized  convection  velocity 
in  his  case  is  0.75  at  an  acoustic  Mach  number  0.9  in  a  cold  jet.  His  results  also  compare  well  with  the 
present  results  at  higher  Stop’s.  In  a  forced  jet  by  glow  discharge,  the  normalized  convection  velocity 
measured  by  Troutt  and  McLaughlin  [1982]  is  about  0.8  for  a  Stop  of  0.3  -0.8  in  a  Mach  2.1  jet.  Although 
their  results  compare  well  with  the  present  values  at  higher  Strouhal  numbers,  they  did  not  see  an  increase 
in  convection  velocity  as  the  Stop  was  increased  from  0.3  to  0.8,  which  potentially  could  indicate  the  lack 
of  actuator  authority  to  force  the  jet  column  instability. 


Fig.  3.1.16  Convection  velocity  for  three  modes. 


3.2  Mach  1.3  Cold  Supersonic  Jet 

As  in  the  results  in  the  Mach  0.9  cold  jet,  the  performance  of  the  plasma  actuators  is  evaluated  by 
PIV  measurements.  The  centerline  Mach  number  decay  and  the  jet  width  (5),  defined  by  the  full  width  at 
half  maximum  (FWHM)  of  the  streamwise  velocity,  is  used  for  overall  performance  evaluation.  The 
spacing  and  convection  velocity  of  generated  large-scale  structures  is  obtained  from  spatial  cross- 
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correlation  and  these  are  used  for  conditional  sampling  of  PIV  images.  The  large-scale  structures  are 
visualized  by  using  conditionally-sampled  Galilean  velocity  fields  (the  coordinate  systems  are  moving 
with  the  convective  velocity  of  large-scale  structures).  From  this  information  about  large-scale  structure, 
the  role  of  the  generated  structures  in  the  jet  development  is  discussed  extensively.  The  jet  flow  fields  are 
measured  at  Mj  =  1.2  (overexpanded),  1.3  (perfectly  expanded),  and  1.4  (underexpanded).  The  results  for 
the  perfectly  expanded  jet  (Mj  =  1.3)  are  very  similar  to  those  in  Mach  0.9  cold  subsonic  jet  as  will  be 
detailed  later. 

3.2.1  Effects  of  Forcing  Strouhal  Number  on  Overall  Jet  Mixing 

The  results  for  Mach  0.9  subsonic  [Kim  et  al.  2007,  2009]  and  perfectly-expanded  Mach  1.3 
supersonic  [Samimy  et  al.  2007]  jets  showed  that  the  forcing  is  most  effective  at  m  =  ±1.  Thus,  the 
results  at  m  =  ±1  are  used  for  the  evaluation  of  the  effects  of  Stop  numbers  on  the  jet  spreading  for  all 
Mj’s.  Average  streamwise  velocity  contours  for  m  =  ±1  are  shown  in  Fig.  3.2.1  for  three  fully-expanded 
jet  Mach  numbers  of  1.2  (over-expanded),  1.3  (perfectly-expanded),  and  1.4  (under-expanded).  The 
streamwise  velocity  is  scaled  from  -40  m/s  to  the  maximum  for  each  Mach  number.  The  maximum  jet 
velocity  is  about  360,  380,  and  420  m/s  for  Mj  =  1.2,  1.3,  and  1.4,  respectively.  Thus,  no  information  can 
be  gained  from  a  one-to-one  comparison  of  the  colors  in  plots  of  differing  Mach  number. 


Figure  3.2.1  Average  streamwise  velocity  contours  for  various  Stop  numbers  at  three  jet  Mach 
numbers.  The  maximum  velocity  of  the  jet  is  about  360,  380,  and  420  m/s  for  Mach  =  1.2, 1.3,  and  1.4 

jets,  respectively. 


The  baseline/unforced  jets  show  that  the  jet  spreading  is  increased  at  the  off-design  jets  of  Mj  = 
1.2  and  1.4.  The  enhanced  spreading  is  due  to  the  feedback  mechanism  sustained  by  upstream-traveling 
acoustic  waves  and  downstream  traveling  large-scale  structures/hydrodynamic  waves  in  the  jet  shear 
layers  that  interact  with  the  shock  waves  generating  the  acoustic  waves.  A  strong  tone  is  also  generated  by 
the  feedback  mechanism  in  all  three  cases,  as  shown  in  Figure  3.2.2.  For  the  imperfectly-expanded  jets. 
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the  broadband  shock  associated  noise  (broad  humps  in  the  spectra)  is  significantly  increased.  However, 
the  shock  associated  broad  noise  is  not  significant  in  the  perfectly  expanded  jet  (Fig.  3.2.2b).  The  shock 
cell  patterns  are  clearly  seen  in  the  average  streamwise  velocity  contours  for  the  imperfectly-expanded  jet 
(Figures  3.2.1a  &  i).  Thus,  the  shock  strength  is  less  in  the  perfectly-expanded  jet  than  in  the  imperfectly- 
expanded  jet.  This  can  be  more  clearly  observed  in  the  centerline  Mach  number  to  be  presented  later. 


(a)  Mj=1.2 


(c)  Mj=1.4 


Fig.  3.2.2  Average  spectra  at  Mj  =  1.2, 1.3  and  1.4,  measured  at  90°  relative  to  the  jet  centerline. 


For  the  over-expanded  jet  (Mj  =  1.2),  the  effect  of  forcing  is  not  apparent  at  a  low  SIdf  of  0.13. 
The  maximum  spreading  occurs  at  a  Stop  of  0.33  (Fig,  3.2.1c),  but  the  enhancement  of  jet  spreading  is 
moderate.  At  a  higher  Stop  of  1.3  (Figure  3.2. Id),  it  appears  that  the  jet  spreading  is  even  suppressed.  The 
contours  for  the  under-expanded  jet  (Figures  3.2.1  i-1)  show  that  the  trend  of  jet  spreading  with  Stop  is 
roughly  similar  to  that  for  the  over-expanded  jet.  For  this  flow  regime,  the  maximum  spreading  is  at  a 
slightly  low  Stop  of  0.27  (Figure  3.2.1k).  As  will  be  further  discussed,  the  forcing  is  less  effective  in  the 
imperfectly-expanded  jets  when  compared  to  the  perfectly-expanded  Mach  1.3  jet. 
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(a)  Mj  =  0.9  [Kim  et  al.  2009]  (b)  Mj  =  1.3 

Figure  3.2.3  Average  velocity  contours  in  the  flapping  plane  at  m  =  ±1  and  Stop  «  0.3.  The  scale  is 
about  the  same,  but  the  spans  in  streamwise  and  cross-streamwise  directions  are  different. 

For  the  perfectly  expanded  jet  (Mj  =  1.3),  the  jet  responds  to  the  forcing  in  a  wide  range  of  StDp‘s. 
At  a  low  Stop  of  0.13,  the  jet  spreading  is  significantly  enhanced,  contrary  to  the  imperfectly-expanded 
cases.  The  maximum  spreading  is  observed  at  a  Stop  of  0.33  (Figure  3.2. Ig)  and  the  enhancement  in  the 
jet  spreading  is  dramatic.  At  a  high  Stop  of  1.3,  the  velocity  contour  is  very  similar  to  that  of  the  baseline. 


32 


implying  that  forcing  is  not  effective  at  high  Stores.  The  trend  observed  at  the  perfectly-expanded  Mach 
1 .3  supersonic  jet  is  very  similar  to  what  was  observed  in  a  subsonic  Mj  0.9  jet  [Kim  et  al.  2009a  &  b]. 


Figure  3.2.3  shows  the  streamwise  velocity  contours  measured  by  the  PIV  system  for  Mach  0.9 
and  1.3  jets,  respectively,  at  a  Stop  of  about  0.3  and  at  m  =  ±1.  Note  that  the  color  map  is  not  the  same  - 
the  same  color  does  not  represent  the  same  velocity.  The  jet  exit  velocity  is  about  280  and  380  m/s  for  Mj 
=  0.9  and  1.3,  respectively.  In  both  jets,  the  actuators  have  control  authority  and  the  enhancement  of 
mixing/spreading  (spreading  from  here  on)  is  about  the  same.  As  will  be  further  discussed  in  a  later 
section,  the  nature  and  role  of  generated  structures  in  the  jet  development  are  also  about  the  same. 


x/D  x/D 


(a)  Jet  width  at  low  SIdf^s 


x/D 


(b)  Jet  width  at  high  Stores 


x/D 


(c)  Equivalent  jet  width  at  low  SIdf^s 


(d)  Equivalent  jet  width  at  high 


Fig.  3.2.4  Jet  width  development  at  m  =  ±1  for  various  Stop  numbers. 


The  effects  of  forcing  Strouhal  number  will  be  more  extensively  presented  by  examining  jet 
width  and  jet  centerline  Mach  numbers.  For  Mj  =  1.3,  the  jet  width  development  at  m  =  ±1  is  shown  in 
Fig.  3.2.4  at  various  Stop’s.  The  jet  width  in  Figs.  3.2.4a&b  is  on  the  flapping  plane,  which  shows  the 
effects  of  forcing  Strouhal  number.  The  jet  width  on  the  non-flapping  plane  (not  shown  here)  does  not 
show  any  significant  spreading.  Thus,  the  cross-section  of  the  jet  plume  is  elliptic  at  this  forcing  mode. 
An  equivalent  jet  width,  defined  as  the  geometric  average  of  the  jet  width  in  the  flapping  and  non¬ 
flapping  planes  (square  root  of  the  multiplication  of  two  jet  widths),  is  shown  in  Figs.  3.2.4c&d.  One 
could  then  compare  this  jet  width  with  those  of  other  modes,  which  are  axisymmetric  in  the  average  sense. 
As  in  the  Mj  =  0.9  subsonic  jet  [Kim  et  al.  2009],  the  jet  plume  spreading  was  significantly  enhanced  by 
forcing.  As  the  Stop  number  is  increased,  the  spreading  also  increases  up  to  Stop  -  0.3  as  shown  in  Figs. 
3.2.4a  &  c.  When  the  Stop  number  is  further  increased,  the  jet  spreading  is  decreased  as  shown  in  Figs. 
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3.2.4b  &  d  and  more  visually  in  Fig.  3.2. Ih.  Thus,  the  enhaneement  of  the  jet  width  is  greatest  at  StoF^ 
0.3.  At  high  StDF^s  greater  than  1.31,  the  jet  width  development  is  about  the  same  as  that  of  the  baseline  as 
ean  be  seen  also  in  Figs.  3.2.1e&h.  The  trend  of  jet  spreading  with  Stj^/s  is  more  readily  seen  in  Fig. 
3.2.5,  showing  the  jet  widths  at  x/D  =  10  for  m  =  ±1.  In  the  figure,  the  jet  width  for  the  foreed  eases  is 
normalized  by  that  for  the  baseline.  The  normalized  jet  width  inereases  rapidly  at  approaehing  0.33. 
For  StDF^s  greater  than  0.33,  the  normalized  width  deereased  with  inereasing  StoF^s  as  was  seen  in  Figs. 
3.2.1  and  3.2.4.  The  jet  width  at  Mj  =  0.9  for  m  =  ±1  (Fig.  3.1.4)  is  very  similar  to  that  at  Mj  =  1.3  shown 
in  Fig.  3.2.4.  These  results  show  that  the  performanee  of  the  aetuators  is  about  the  same  in  both  Maeh  0.9 
subsonie  and  perfeetly-expanded  Maeh  1.3  jets. 


Fig.  3.2.5  Normalized  jet  widths  on  the  flapping  plane  at  x/D  =  10  for  m  =  ±1.  The  jet  width  at  each 
StoF  number  was  normalized  by  that  for  the  baseline. 


(a)  Mj=1.2 


0.20 

-0.23 


xfO 

(b)  Mj=1.4 


(c)  Normalized  jet  width  at  x/D  =10 


Fig.  3.2.6  Jet  width  development  on  the  flapping  plane  at  m  =  ±1  for  the  imperfectly-expanded  jets. 


For  the  imperfectly-expanded  Mach  1.2  &  1.4  jets,  the  jet  width  development  with  downstream 
location  at  m  =  ±1  is  shown  in  Figure  3.2.6.  In  both  jets,  the  jet  width  increases  at  low  StoF^s  less  than 
about  0.3.  The  enhancement  of  jet  spreading  is  maximum  at  StoF  numbers  0.33  and  0.26  for  Mj  =1.2  and 
1.4,  respectively.  The  jet  width  shows  a  dip  and  secondary  peak  at  high  StoF  numbers  when  the  StoF 
number  is  increased  further  from  the  maximum  (Fig.  3.2.6c),  which  was  not  seen  in  the  perfectly- 
expanded  jet  (Fig.  3.2.5).  For  some  other  azimuthal  modes  (not  shown  here),  the  normalized  jet  width  is 
undulating  with  StoF  numbers.  This  difference  in  jet  width  trend  is  possibly  due  to  the  interaction  of  the 
forced  and  naturally  amplified  (by  the  feedback  loop)  structures  as  will  be  further  discussed  later.  At 
StDF’s  greater  than  1.0,  the  jet  width  is  reduced  by  forcing  as  was  also  observed  in  the  velocity  contours  in 


34 


Fig.  3.2.1.  The  overall  enhaneement  of  jet  spreading  is  not  as  signifieant  as  in  the  perfeetly-expanded  jets. 
It  seems  that  the  reduetion  in  jet  spreading  at  high  Stop  numbers  and  overall  spreading  is  also  assoeiated 
with  the  interaetion  of  the  foreed  and  naturally  oeeurring  struetures. 


Table  3.  Optimal  Stop  numbers,  showing  maximum  jet  spreading  at  each  azimuthal  mode. 


Azimuthal  mode 

Mj=1.2 

Mj=1.3 

Mj=  1.4 

m  =  0 

0.52 

0.52 

0.33 

m=  1 

0.26 

0.39 

0.26 

m  =  2 

0.26 

0.33 

0.06 

m  =  3 

0.20 

0.13 

0.26 

m  =  ±l 

0.33 

0.33 

0.26 

(N 

-H 

II 

a 

0.52 

0.33 

0.46 

For  other  azimuthal  modes,  the  optimal  SIdf  numbers  are  seleeted  from  the  normalized  jet  width  at 
x/D  =  10  and  are  shown  in  Table  3.  For  the  most  eases,  the  optimal  Stop  number  is  about  0.3. 
Exeeptionally  low  numbers  are  seen  at  0.13  for  Mj  =  1.3  &  m  =  3  and  0.06  for  Mj  =  1.4  &  m  =  2.  In  the 
Mj  =  0.9  subsonie  jet,  the  optimal  SIdp  number  was  lower  than  any  other  modes  as  in  the  Mj  =  1.3 
perfeetly-expanded  supersonie  jet.  For  other  eases,  the  numbers  are  within  0.2-0. 6  range,  found  in  the 
literature. 

3.2.2  Effects  of  azimuthal  modes 

3.2.2. 1  Perfectly-expanded  jet 

The  results  presented  in  the  earlier  seetion  showed  the  effeets  of  Stop  numbers  at  m  =  ±1.  In  this 
seetion,  the  effeets  of  azimuthal  modes  will  be  diseussed  by  using  optimal  eases;  those  that  show  the  most 
jet  spreading.  For  the  perfeetly-expanded  jet,  the  average  streamwise  veloeity  eontours  at  the  optimal  SIdf 
numbers,  listed  on  Table  2,  are  shown  in  Fig.  3.2.7.  Also,  the  profiles  of  the  eenterline  Maeh  number  and 
jet  width  are  shown  in  Figure  3.2.8  for  the  optimal  SIdf^s.  Note  that  an  equivalent  jet  width  is  used  only 
for  m  ±1  sinee  the  jet  eross-seetion  is  elliptie  for  this  mode.  The  streamwise  veloeity  eontours  and 
eenterline  Maeh  number  show  that  the  jet  potential  eore  length  is  shortened  signifieantly  for  m  =  1  and  ±1, 
and  moderately  for  the  other  modes.  The  potential  eore  length  is  redueed  from  7  nozzle  diameters  in  the 
baseline  jet  to  4  (for  m  =  1  &  ±1)  and  5.5  (for  the  rest  of  modes).  Although  both  the  potential  eore  length 
and  eenterline  Maeh  number  are  indireet  measures  for  the  jet  growth/spreading,  the  trend  observed  in  Figs. 
3.2.7  and  3.2.8a  is  very  similar  to  what  is  seen  in  the  jet  width,  a  direet  measure  for  the  spreading.  There 
is  a  moderate  undulation  in  the  eenterline  Maeh  number  due  to  weak  shoek  eells.  Although  the  diverging 
seetion  of  the  nozzle  was  designed  by  the  method  of  eharaeteristies,  the  oeeurrenee  of  weak  shoeks  is 
unavoidable  with  a  thiek  lipped  nozzle.  For  all  azimuthal  modes,  the  eenterline  Maeh  number  deeay  is 
enhaneed  by  foreing.  The  results  in  Figs.  3.2. 7-8  show  that  the  most  effeetive  foreing  is  at  m  =  ±1  as  in 
the  Mj  =  0.9  subsonie  jet  (See.  3.1).  At  this  mode,  a  dramatie  enhaneement  in  jet  spreading  is  manifested 
in  aeeelerated  eenterline  Maeh  number  deeay  and  enhaneed  jet  width.  Additionally,  the  spreading  at  m  = 
1  is  substantially  improved.  For  the  rest  of  the  azimuthal  modes,  the  enhaneement  in  jet  spreading  is 
moderate. 
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Fig.  3.2.7  Average  streamwise  velocity  contours  at  the  optimal  Stores,  listed  in  Table  2,  for  each 

azimuthal  mode. 


(a)  Centerline  Maeh  number 

Fig.  3.2.8  Comparison  of  the  centerline  Mach  number  and  jet  width  at  Mj  =  1.3.  For  m  =  ±1,  the  jet 
width  is  the  equivalent  width.  The  forcing  Strouhal  number  is  shown  in  the  legend  of  each  figure. 


The  results  in  the  perfeetly-expanded  Mj  1.3  jet  are  largely  similar  to  those  in  Mj  0.9  subsonie  jets 
[Kim  et  al.  2009].  However,  the  jet  growth  at  m  =  ±2  is  signifieantly  redueed  eompared  to  those  in  Maeh 
0.9  jets.  In  the  Mj  0.9  jet,  the  inerease  in  eenterline  Maeh  number  deeay  was  similar  for  both  m  =  ±2  and 
m  =  ±1  eases.  In  the  perfeetly-expanded  Maeh  1.3  jet,  the  jet  growth  at  m  =  ±2  is  about  the  same  as  the 
moderately  effeetive  group  of  modes  (m  =  0,  &  2-3).  At  this  point,  it  is  not  elear  why  the  performanee  at 
m  =  ±2  was  redueed  in  Mj  1.3  jets.  Another  differenee  is  the  growth  in  the  initial  shear  layer.  In  the  Mj 
0.9  subsonie  jet,  a  signifieant  enhaneement  in  jet  width  in  the  initial  shear  layer  was  seen  for  many 
azimuthal  modes  at  a  relatively  high  Stop  of  about  1.0  (see  Fig.  3.1.7)  [Kim  et  al.  2007,  2009].  For  the 
perfeetly-expanded  Mj  1.3  jet,  the  enhaneement  in  jet  width  in  the  initial  shear  layer  is  only  seen  at  m  =  3 
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(shown  in  Fig.  3.2.9)  and  is  not  as  significant  as  in  the  subsonic  counterpart.  The  increased  jet  width 
upstream  of  the  end  of  potential  core  is  thought  to  be  due  to  growths  of  the  structures  generated  by 
forcing. 


x/D 

Fig.  3.2.9  Jet  width  development  along  the  streamwise  direction  at  m  =  3. 


(a)  Centerline  Mach  number  (b)  Jet  width 

Fig.  3.2.10  Comparison  of  the  centerline  Mach  number  and  jet  width  at  Mj  =  1.2.  The  jet  width  for 
m  =  ±1  is  on  the  flapping  mode.  The  forcing  Strouhal  number  is  0.52,  0.26,  and  0.33  for  m  =  0, 1, 

and  ±1,  respectively. 


3.2.2.2  Imperfectly-expanded  jets 

The  centerline  Mach  number  and  jet  width  development  for  two  modes  of  m  =  0  and  1  are 
compared  to  show  the  effects  of  azimuthal  mode  in  the  over-expanded  jets  (Fig.  3.2.10).  In  this  figure,  the 
Stop  number  for  each  azimuthal  mode  is  selected  for  the  maximum  jet  spreading  as  listed  in  Table  3.  In 
the  over-expanded  Mj  =  1.2  jet,  the  centerline  Mach  number  is  a  better  measure  for  the  overall  jet 
spreading  since  no  measurement  is  done  on  the  non-flapping  plane  for  the  m  =  ±1  case.  Neither  the 
potential  core  length,  obtained  from  Fig.  3.2.10a,  nor  the  centerline  Mach  number  decay  rate  beyond  the 
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end  of  potential  eore  is  signifieantly  altered  by  foreing.  The  ineffeetiveness  of  foreing  in  the  over- 
expanded  jet  is  also  manifested  in  the  average  streamwise  veloeity  eontours  shown  in  Fig.  3.2.11.  The 
shoek  eell  patterns  in  the  potential  eore  region  and  jet  growth  are  barely  ehanged  by  foreing.  However, 
the  jet  width  on  the  flapping  plane  at  m  =  ±1  shows  a  notable  inerease,  suggesting  that  large-seale 
struetures  are  generated.  However,  it  seems  that  the  generated  large-seale  struetures  laek  the  strength  to 
signifieantly  exeite  the  shoek  eontaining  jet  and  so  do  not  inerease  mixing.  This  will  to  be  diseussed 
further  later.  For  the  other  modes  not  presented  here,  the  spreading  is  about  the  same  as  that  for  m  =  0  or 
1.  These  results  show  that  the  foreing  in  the  over-expanded  jet  is  not  as  effeetive  as  in  the  perfeetly- 
expanded  jet. 


For  the  under-expanded  Mj  =  1.4  jet,  the  eenterline  Maeh  number  and  jet  width  are  shown  in  Fig. 
3.2.12  for  m  =  0,  1,  and  ±1.  For  the  other  modes  not  shown  here,  the  eenterline  Maeh  number  and  jet 
width  are  very  similar  to  that  for  m  =  1  and  m  =  0,  respeetively.  Again,  in  this  figure,  the  foreing  Strouhal 
number  for  eaeh  azimuthal  mode  is  seleeted  for  the  maximum  jet  spreading  as  listed  in  Table  3.  The 
eenterline  Maeh  number  undulates  between  1.3  and  1.5  due  to  periodie  shoek  eell  struetures  in  the  jet. 
The  jet  potential  eore  length  and  eenterline  Maeh  number  deeay  are  not  signifieantly  ehanged  by  foreing, 
as  in  the  over-expanded  jet.  The  jet  width  is  inereased  slightly  by  foreing,  but  the  inerease  is  not  as 
signifieant  as  in  Mj  =  0.9  subsonie  (See.  3.1)  or  ideally  expanded  Mj  =  1.3  jets.  The  jet  width 
enhaneement  at  m  =  ±1  seems  signifieant,  but  note  that  the  non-flapping  plane  width,  whieh  would  be 
very  elose  to  that  for  the  baseline,  is  not  taken  into  aeeount.  The  equivalent  jet  width  is  expeeted  to  be 
very  elose  to  that  for  the  m  =  0.  This  explains  why  the  eenterline  Maeh  number  deeays  for  m  =  1  and  ±1 
are  very  elose  to  eaeh  other  as  seen  in  Fig.  3.2.12a.  The  redueed  jet  spreading  over  the  baseline  jet  is 
partially  due  to  enhaneed  mixing  in  the  baseline/unforeed  jet  as  diseussed  earlier  (Figs.  3.2.1a,e,&i). 

In  Mj  =  1.4  jets,  the  eenterline  Maeh  number  deeay  is  slightly  suppressed  for  all  azimuthal  modes 
exeept  for  m  =  ±1,  as  shown  in  Fig.  3.2.12a.  For  m  =  0,  the  eenterline  Maeh  number  deeay  is  signifieantly 
redueed  and  the  jet  width  development  is  almost  the  same  as  the  baseline.  In  the  Mj  0.9  subsonie  jet  (See. 
3.1),  it  was  shown  that  vortex  rings  were  generated  at  m  =  0,  and  that  the  eenterline  veloeity  deeay  and  jet 
spreading  were  redueed  due  to  self  induetion  and  axisymmetrie  nature  of  the  vortex  ring  [Kim  et  al.  2007 
&2009]. 
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(a)  Centerline  Maeh  number  (b)  Jet  width 

Fig.  3.2.12  Comparison  of  centerline  Mach  number  and  jet  width  at  Mj  =  1.4.  The  jet  width  for  m 
=  ±1  is  on  the  flapping  mode.  The  number  in  the  parentheses  indicates  the  forcing  Strouhal 

number  for  each  azimuthal  mode. 


3.2.3  Effects  of  Forcing  on  Turbulence 

This  section  examines  the  development  of  jet  centerline  two-dimensional  turbulent  kinetic  energy 
(TKE  henceforth),  as  only  a  two-component  PfV  system  was  utilized.  The  effects  of  forcing  Strouhal 
number  on  TKE  is  shown  in  Fig.  3.2.13  in  the  perfectly-expanded  jet.  The  TKE  level  is  significantly 
increased  at  a  Stop  number  near  of  0.3,  but  its  level  is  close  to  that  for  the  baseline  at  low  (not  shown  here) 
and  high  Stop  numbers  (1.31,  for  an  example).  For  the  Stop  of  0.33  which  showed  maximum  jet  growth, 
the  TKE  level  saturates  at  x/D  =  8  and  then  slowly  decays.  At  other  Stop’s,  the  TKE  level  increases 
almost  monotonically  without  showing  any  saturation  in  the  entire  streamwise  measurement  span.  This 
trend  was  not  observed  in  the  MJ  0.9  subsonic  jet,  where  TKE  saturated  at  x/D  =  7  at  low  Stop’s  (Fig. 
3.1.10). 


Fig.  3.2.13  Turbulent  kinetic  energy  development  at  m  =  ±1  in  the  perfectly-expanded  jet. 


Figure  3.2.14  shows  the  effects  of  azimuthal  modes  for  three  flow  regimes  of  Mj  =  1.2,  1.3,  and 
1.4.  The  cases  shown  in  the  figure  are  for  forcing  Strouhal  numbers  that  achieved  maximum  spreading  for 
each  mode,  listed  in  Table  3.  The  forcing  Strouhal  numbers  for  each  mode  correspond  to  those  shown 
Figs.  3.2.8,  3.2.10,  and  3.2.12,  respectively.  For  Mj  =  1.3,  TKE  is  significantly  increased  for  all  forcing 
modes  as  in  the  Mj  0.9  subsonic  jet.  The  increase  in  TKE  level  in  the  Mj  1.3  supersonic  jet  is  a  little  bit 
less  than  that  in  the  Mj  0.9  subsonic  jet  (compare  Figs.  3.1.10a  and  3.2.14a).  The  TKE  saturates  at  x/D  = 
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Normalized  TKE 


8.5  and  10.5  for  m  =  ±1  and  1,  respectively.  This  suggests  that  the  earlier  saturation  of  TKE  indicates 
increased  jet  growth/spreading.  As  in  a  Mj  =  0.9  subsonic  jet,  the  TKE  is  still  on  the  rise  for  m  =  0.  This  is 
due  to  self-induction  by  vortex  rings  and  the  symmetric  nature  of  the  generated  structures  in  this  mode 
[Kim  et  al.  2007  &  2009],  as  discussed  in  Sec.  3.1. 


Fig.  3.2.14  Two-dimensional  turbulent  kinetic  energy  development  for  three  azimuthal 

modes  m  =  0, 1,  and  ±1. 

For  imperfectly-expanded  jets  of  Mj  =  1.2  and  1.4,  it  appears  that  centerline  TKE  is  not 
significantly  altered  by  forcing  as  seen  in  Fig.  3.2.14b  &c.  This  is  partially  due  to  increased  TKE  level  in 
the  unforced  baseline  jets  due  to  the  presence  of  relatively  strong  shock/expansion  train.  As  was  discussed 
earlier,  when  the  jet  operates  in  imperfectly-expanded  regime,  a  screech  tone  is  generated  by  a  feedback 
mechanism  between  the  hydrodynamic  and  acoustic  waves  (Fig.  3.2.2),  and  large-scale  structures  are 
amplified  due  to  this  feedback  mechanism  resulting  in  the  increase  in  TKE.  If  the  large-scale  structures 
are  suppressed,  TKE  would  be  reduced  by  forcing.  This  was  observed  when  the  jet  was  forced  at  high 
StDF^s  (not  shown  here). 

3.2.4  Large-Scale  Structures  and  their  role  in  the  Jet  Development 

Large-scale  structures  are  visualized  by  using  conditionally  sampled  Galilean  velocity  field.  In 
the  Galilean  velocity  field,  the  reference  frame  moves  with  the  convection  velocity  of  large-scale 
structures  in  the  flow.  Thus,  the  large-scale  structures  are  stationary  in  this  frame,  and  they  are  identified 
if  the  streamlines  show  closed  or  spiral  shapes  [Kline  and  Robinson  1990,  Robinson  et  al.  1989].  The 
procedure  was  discussed  in  more  detail  earlier. 

3.2.4.1  Perfectly-Expanded  Jet  (Mj  =  1.3) 

Large-scale  structures  for  excitation  with  m  =  ±1  at  various  Strouhal  numbers  in  the  Mj  =  1.3  jet 
are  shown  Fig.  3.2.15.  The  relative  magnitude  of  the  streamwise  velocity  is  represented  by  color:  red  and 
blue  indicate  fast  and  slow  speeds,  respectively.  There  is  no  common  color  map,  valid  for  all  images, 
since  the  reference  frame  moves  at  a  different  speed  for  each  case.  The  structures  identified  by  Galilean 
streamlines  are  similar  to  those  visualized  by  laser  scattering  [Samimy  et  al.  2007b].  Large  periodic 
structures  are  generated  by  forcing  at  a  wide  range  of  St^/s  from  0.2  to  1.05  as  in  the  Mj  0.9  subsonic  jet 
(Sec.  3.1).  At  StDF^s  outside  of  this  range,  there  are  no  visible  periodic  structures  in  the  shear  layer  as 
shown  in  Fig.  3.2. 15e,  as  an  example.  When  this  figure  is  compared  to  the  case  for  Mj  0.9  jet  (Fig.  3.1.11), 
the  jet  spreading  and  structure  formation  by  forcing  are  very  similar  in  qualitative  sense. 
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(d)  Stop  =  1 .05  (e)  Stop  =1.31  (f)  Baseline 

Fig.  3.2.15  Galilean  streamlines  superimposed  on  the  streamwise  velocity  fields  at  m  =  ±1  in  the  Mj  = 
1.3  jet.  The  velocity  fields  were  conditionally-averaged  and  the  number  of  images  used  for  the 

averaging  was  30. 

At  SIdp  =  0.33,  the  generated  stmetures  are  very  robust  and  well  organized.  The  generated 
struetures  are  nearly  eireular  in  shape,  penetrating  into  the  jet  eenterline,  and  eausing  signifieant 
undulation  in  the  jet  plume.  The  potential  eore  length  is  about  four  nozzle  exit  diameters  as  shown  in  Fig. 
3.2.8a.  The  potential  eore  is  signifieantly  shortened  by  the  entrainment  of  ambient  air  and  penetration  of 
energetie  struetures  into  the  jet  as  seen  in  Fig.  3.2.15b.  Also,  the  jet  spreading  is  signifieantly  enhaneed  by 
the  robust  struetures  as  shown  in  Fig.  3.2.8b.  The  generated  struetures  either  deeays  eompletely  and/or 
beeome  very  disorganized  by  x/D  =  7  -  no  identifiable  large  struetures  are  seen  downstream  of  this 
loeation.  As  shown  in  Fig.  3.2.13,  the  TKE  development  is  saturated  at  this  loeation  and  this  saturation 
loeation  is  related  to  the  sudden  deeay  of  large-seale  struetures. 

At  a  low  Stop  number  of  0.2,  the  generated  large  periodie  struetures  appear  to  be  nearly  elliptieal 
with  a  high  aspeet  ratio  aligned  in  the  streamwise  direetion.  This  may  be  due  to  a  limited  growth  of  the 
struetures  in  the  spanwise  direetion  when  eompared  to  that  at  Stop  =  0.33.  Also,  the  interaetion  between 
the  generated  struetures  and  the  jet  plume  appears  to  be  limited.  The  monotonie  inerease  in  TKE  (Fig. 
3.2.13)  up  to  x/D  of  12  indieates  that  the  interaetion  between  struetures  is  signifieantly  less  for  this  ease. 
This  ean  be  inferred  from  a  longer  spatial  lifetime  of  the  generated  struetures  at  Stop  =  0.2  than  at  Stop  = 
0.33.  Some  large  struetures  are  seen  up  to  x/D  =  10  at  Stop  =  0.2  while  all  struetures  had  deeayed  by  x/D 
=  7  for  StDP  =  0.33.  As  in  Stop  =  0.33  ease,  the  jet  plume  is  undulating  in  the  lateral  direetion  due  to  the 
flapping  aetion  eaused  by  the  generated  struetures.  However,  the  undulating  motion  is  not  as  signifieant 
as  in  Stop  =  0.33  ease  due  to  the  less  energetie  struetures.  As  a  result,  the  mixing  enhaneement  is 
relatively  less  signifieant  at  this  foreing  Strouhal  number. 

At  a  moderate  Stop  of  0.52,  periodie  large  struetures,  with  a  redueed  spaeing,  are  generated.  The 
interaetion  between  generated  struetures  seems  minimal,  whieh  is  inferred  from  the  well-preserved 
periodie  struetures  and  their  almost  eonstant  dimension  up  to  x/D  =  8.  It  seems  that  the  jet  plume 
undulates  less  sinee  the  generated  struetures  are  smaller  than  those  at  low  Stop^s.  The  redueed  undulation 
of  the  jet  plume  suggests  that  the  interaetion  between  the  generated  struetures  and  the  jet  plume  is  not 
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significant  when  compared  to  that  at  low  Stores.  The  reduced  TKE  at  this  Stop  number,  shown  in  Fig. 
3.2.13,  also  confirms  that  the  interaction  was  decreased  or  limited  at  this  moderate  Stop. 


Interestingly,  some  periodic  structures  are  observed  even  at  a  high  Stop  number  of  1.05  (Fig. 
3.2. 15d).  The  generated  structures  are  very  small  and  closely  spaced.  It  appears  that  these  small 
structures  decay  faster  than  those  at  low  Stop  numbers,  but  they  are  occasionally  visible  up  to  x/D  =  8.  At 
a  higher  Stop  number  of  1.3  (Fig.  3.2. 15e),  no  periodic  structures  are  seen  and  the  flow  fields  are  very 
similar  to  that  of  the  baseline  (Fig.  3.2. 15f).  In  the  baseline  jet,  there  are  some  randomly  spaced  structures, 
but  they  are  not  energetic  so  that  their  effect  on  the  jet  plume  is  minimal.  The  structures  observed  in  this 
Mj  =  1.3  jet  are  very  similar  to  what  were  seen  in  the  Mj  =  0.9  subsonic  jet  with  a  Reynolds  number  of 
about  0.7  X  10^  (Sec.  3.1). 


Figure  3.2.16  shows  the  effects  of  Stop  number  on  the  spacing  (and  also  dimension)  of  the 
generated  structures.  The  structure  spacing  is  calculated  from  two-dimensional  spatial  correlation  of  700 
instantaneous  velocity  fields  as  detailed  in  Kim  et  al.  [2009a&b].  The  spacing  (X.)  is  inversely 
proportional  to  Stop  number  as  shown  in  Fig.  3.2.16,  and  the  streamwise  dimension  of  the  structures  also 
shows  the  same  trend  as  the  structure  spacing,  as  seen  in  Fig.  3.2.15.  The  profiles  for  various  azimuthal 
modes  in  Fig.  3.2.16  are  collapsed  into  a  single  curve  using  the  following  equation 


X 

D 


Stop 


+  C 


(3) 


where  a  and  c  are  constants.  The  figure  shows  that  the  structure  spacing  is  strongly  dependant  on  the  Stop 
number,  and  that  effects  of  azimuthal  modes  are  minimal.  Also  shown  in  the  figure  is  that  the  perfectly- 
expanded  jet  responds  to  the  forcing  over  a  range  of  Stop  numbers  from  0.2  to  1.3.  This  strong  relation 
further  confirms  that  the  structures  seen  in  Fig.  3.2.15  were  generated  by  the  excitation/forcing  rather  than 
by  any  other  means. 


Fig.  3.2.16  Spacing  of  the  generated  structures  in  Mj  =  1.3  jets. 


Figure  3.2.17  shows  visualized  large-scale  structures  for  m  =  0  and  1  at  Stop  numbers  showing  the 
maximum  spreading  for  each  mode  (Fig.  3.2.8  and  Table  3).  For  m  =  0,  the  visualized  structures  are 


42 


symmetric  across  the  jet  plume  since  they  are  vortex  rings.  The  centerline  streamwise  velocity  between  a 
pair  of  vortical  structures  (actually  inside  a  vortex  ring)  is  faster  due  to  self  induction,  but  the  flow 
between  two  neighboring  vortex  rings  is  slower  due  to  the  entrainment  of  slow  moving  ambient  air  [Kim 
et  al.  2009a&b].  Thus,  the  streamwise  velocity  along  the  jet  centerline  undulates,  similar  to  what  was 
observed  in  Mj  =  0.9  subsonic  jet  [Kim  et  al.  2009a].  As  in  the  Mach  0.9  jet,  the  structures  do  not 
penetrate  into  the  jet  centerline  due  to  the  symmetric  nature  of  the  generated  structures.  Consequently, 
the  interaction  of  large-scale  structures  and  the  jet  plume  is  less  destructive.  As  a  result,  the  centerline 
Mach  number  decay  is  less  than  in  that  of  any  other  mode  as  shown  in  Fig.  3.2.8a.  For  m  =  1,  the  pattern 
of  the  generated  structures  and  jet  plume  undulation  is  very  similar  to  that  of  m  =  ±  1,  but  the  structures 
seem  less  energetic.  As  a  result,  the  enhancement  in  the  jet  width  at  this  mode  is  less  than  that  of  m  =  ±1 
(Fig.  3.2.8b). 


Fig.  3.2.17  Galilean  streamlines  superimposed  on  the  conditionally-averaged  streamwise  velocity 

fields  for  the  cases  shown  in  Fig.  3.2.8  or  Table  3. 

For  other  modes  not  shown  here,  the  generated  vortex  pattern  is  very  similar  to  that  for  m  =  0  and 
1  for  even-  and  odd-numbered  modes,  respectively.  For  example,  the  vortex  pattern  of  m  =  2  is  similar  to 
that  at  m  =  0,  shown  in  Fig.  3.2.17a.  In  Mach  0.9  subsonic  jet,  the  streamwise  component  of  turbulence 
was  more  amplified  at  even-numbers  modes,  while  the  cross-streamwise  component  was  more  amplified 
at  odd-numbered  modes.  As  discussed  in  Sec.  3.1,  this  is  due  to  symmetric  or  asymmetric  nature  of  the 
generated  structures  for  even-  or  odd-numbered  modes,  respectively  [Kim  et  al.  2009a&b].  All  these 
results  show  that  the  effects  of  forcing  on  the  flow  structures  and  jet  development  are  very  similar  to  what 
was  observed  in  Mj  =  0.9  subsonic  jet  [Kim  et  al.  2009a].  Thus  the  discussion  of  the  role  of  generated 
structures  on  the  jet  development  presented  in  the  subsonic  case  is  still  valid  in  this  perfectly-expanded 
Mj=  1.3  supersonic  jet. 


3.2.4.2  Imperfectly-Expanded  Jets  (Mj  =  1.2  and  1.4) 

Flow  visualizations  were  conducted  at  Mj  =  1.2  and  1.4  to  investigate  whether  the  plasma 
actuators  are  effective  in  forcing  flows  containing  a  shock/expansion  train.  Some  preliminary  results 
presented  in  Samimy  et  al.  [2007b,  2008]  and  the  results  discussed  in  earlier  sections  showed  that  the 
forcing  is  less  effective  in  imperfectly-expanded  jets  (Figs.  3.2.8,  3.2.10,  &  3.2.12).  The  flow 
visualizations  based  on  the  condensed  water  particle  showed  that  the  jet  responded  to  the  forcing  in  a 
similar  fashion  as  in  perfectly-expanded  Mj  =  1.3  jet  [Samimy  et  al.  2008].  However,  it  seemed  that  the 
jet  did  not  respond  to  the  actuation  in  the  over-expanded  Mj  =  1.2  jet.  The  generated  structures  are 
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visualized  based  using  the  Galilean  veloeity  field  to  find  an  answer  for  the  redueed  effeetiveness  at  the 
imperfeetly-expanded  j  ets. 


Visualized  large-seale  struetures  are  shown  in  Fig.  3.2.18  for  m  =  ±1  in  the  over-expanded  Mj  = 
1.2  jet.  These  struetures  were  also  seen  in  flow  visualizations  in  the  earlier  researeh  [Samimy  et  al.  2008]. 
The  generated  struetures  seem  as  robust  and  energetie  as  in  the  perfeetly-expanded  jet  (Fig.  3.2.15). 
However,  the  enhaneement  in  jet  growth  is  not  as  signifieant  as  that  in  Mj  =  1.3  jets  as  seen  in  Fig. 
3.2.10b.  The  spaeing  of  the  struetures,  the  distanee  between  two  eonseeutive  spiral  shapes,  are  inversely 
proportional  to  Stop  numbers  ranging  from  0.2  to  1.1  as  also  shown  in  Fig.  3.2.19.  For  m  =  0  and  1,  the  jet 
responds  in  a  narrower  range  of  Stop  numbers  from  0.3  to  0.8  when  eompared  to  the  perfeetly-expanded 
jet  ease  shown  in  Fig.  3.2.16.  In  the  baseline  jet,  some  periodie  struetures  are  observed,  but  they  appear  to 
be  not  well  organized.  However,  an  image  aequired  from  the  proper  orthogonal  deeomposition  (not 
shown  here)  showed  that  there  are  periodie  struetures  in  the  baseline  jet.  The  periodie  struetures  in  the 
proper  orthogonal  deeomposition  images  do  not  neeessarily  indieate  that  there  are  sueh  struetures  in  the 
flow  in  steady  fashion.  This  suggests  that  the  spatially  periodie  struetures  are  generated  by  a  feedbaek 
meehanism  in  the  baseline,  but  that  they  are  not  steady  in  time. 


(e)  Stop  =  0.52 


(d)  Baseline 


Fig.  3.2.18  Galilean  streamlines  superimposed  on  the  conditionally-averaged  streamwise  velocity 

fields  at  m  =  ±1  in  the  over-expanded  Mj  =  1.2  jet. 


There  was  a  sign  of  eompetition  for  energy  between  generated  struetures  due  to  foreing  and 
naturally  amplified  struetures  at  a  Stop  about  0.3.  The  indieator  of  this  eompetition  is  the  behavior  of 
strueture  spaeing  espeeially  for  m  =  0.  There  is  dual  spaeing  of  struetures  around  a  Stop  of  0.3  whieh 
suggests  eompetition.  The  presenee  of  dual  spaeing  at  a  Stop  about  0.3  is  observed  in  Fig.  3.2.18b.  In  the 
upstream  region,  the  spaeing  is  very  elose  to  that  for  the  baseline.  However,  the  spaeing  in  the 
downstream  region  eollapses  on  the  relation  in  Eq.  3,  implying  that  the  struetures  are  generated  by  foreing. 
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At  high  Stores  near  0. 8-1.0,  the  foreed  stmetures  are  small  and  less  energetie  so  that  the  naturally 
amplified  stmetures  seem  to  prevail  -  inferred  from  the  stmeture  spaeing.  Unlike  in  the  perfeetly- 
expanded  jet,  the  foreed  stmetures  need  to  eompete  with  the  naturally  amplified  stmetures  to  survive.  As 
a  result  of  this  eompetition  or  interaetion,  it  is  thought  that  the  vortiees  are  weaker  than  in  the  perfeetly- 
expanded  jet  although  it  is  not  readily  seen  in  the  visualized  stmetures  in  Fig.  3.2.18.  A  eareful 
eomparison  of  Figs.  3.2.15  and  3.2.18  suggests  that  the  generated  stmetures  at  this  Stop  are  perhaps  well 
organized,  but  the  spanwise  dimension  is  smaller  than  that  observed  in  the  perfeetly-expanded  Mj  =  1.3 
jet.  These  findings  suggest  that  the  generated  stmeture  in  the  over-expanded  jet  is  less  energetie  than 
those  in  the  perfeetly-expanded  jet.  This  redueed  strength  and  growth  of  the  stmetures  may  be  responsible 
for  the  deereased  effeetiveness  in  mixing  enhaneement  in  the  over-expanded  jet  (Fig.  3.2.10). 

For  the  under-expanded  Mj=1.4  jet,  visualized  stmetures  at  m  =  ±1  are  shown  in  Fig.  3.2.20.  As 
was  observed  in  other  flow  regimes,  the  spaeing  and  dimension  of  stmetures  deereased  with  inereasing 
Stop  number.  In  the  baseline  jet,  there  are  periodie  stmetures  generated  by  a  feedbaek  meehanism  as  in  the 
over-expanded  jet.  The  pattern  of  vortieal  stmetures  in  foreed  eases  is  very  similar  to  what  was  seen  in  the 
other  two  flow  regimes  of  Mj  =  1.2  and  1.3  (Figs.  3.2.15  &  3.2.18).  At  low  SIdp^s  of  0.13  and  0.26,  the 
spaeing  in  the  vieinity  of  the  nozzle  exit  is  smaller  than  that  at  far-downstream  loeations.  The  spaeing  in 
the  upstream  region  is  aetually  very  elose  to  that  for  the  baseline.  As  in  the  over-expanded  jet,  this  dual 
spaeing  in  the  shear  layer  suggests  that  the  foreed  stmetures  need  to  eompete  with  the  naturally  amplified 
stmetures  whieh  ean  be  observed  in  Fig.  3.2.20f  At  moderate  Stop  numbers  of  0.39  and  0.52,  only  single 
spaeing  is  seen  over  the  entire  streamwise  span.  A  dual  spaeing  is  also  observed  at  high  Stop  numbers 
about  1.0,  but  it  is  not  readily  observed  in  Fig.  3.2.20e. 


Fig.  3.2.19  Spacing  of  the  structures  in  Mj  =  1.2  jets. 

Although  the  generated  stmetures  appear  to  be  as  energetie  as  those  in  the  perfeetly-expanded 
ease,  it  is  expeeted  that  the  strength  of  the  generated  stmetures  will  be  less  in  the  under-expanded  jets  due 
to  eompetition.  As  in  the  over-expanded  ease,  the  redueed  strength  may  be  partially  responsible  for  the 
redueed  jet  growth  as  shown  in  Fig.  3.2.12b.  At  Stop  =  0.26,  the  jet  responded  to  the  foreing  and  the 
generated  stmetures  seem  to  be  more  energetie  eompared  to  other  eases  shown  in  the  figure.  These 
energetie  stmetures  are  responsible  for  better  mixing  enhaneement  as  shown  in  Fig.  3.2.12b.  At  Stop  = 
0.39  and  0.52,  the  spaeing  of  the  generated  stmetures  is  very  elose  to  that  of  the  shoek  eells.  The 
generated  stmetures  are  well  organized,  but  it  seems  that  they  are  not  suffieiently  energetie  to  be  able  to 
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undulate  the  jet  eolumn  signifieantly,  as  inferred  from  the  relatively  straight  jet  plume.  When  the  jet  is 
foreed  at  high  as  in  Fig.  3.2.20e,  the  strueture  spaeing  is  the  same  as  the  baseline,  but  the  formation 
of  struetures  is  suppressed  by  foreing  as  inferred  from  the  dimension  and  irregular  spaeing  of  the 
struetures.  The  suppression  of  naturally  amplified  struetures  is  most  likely  responsible  for  the  redueed  jet 
growth  and  the  inereased  potential  eore  length  at  high  Stores  (Fig.  3.2.12b). 


(d)  StDF  =  0.52  (e)  StDF=1.31  (f)  Baseline 

Fig.  3.2.20  Galilean  streamlines  superimposed  on  the  conditionally-averaged  streamwise  velocity 
fields  at  m  =  ±1  in  the  under-expanded  Mj  =  1.4  jet. 


Figure  3.2.21  shows  the  spaeing  of  periodie  struetures  either  foreed  or  naturally  amplified.  The 
spaeing  of  generated  struetures  is  inversely  proportional  to  Stop  numbers  ranging  roughly  from  0.3  to  1.3. 
For  low  and  high  SIdf  numbers,  the  strueture  spaeing  is  the  same  as  that  in  the  baseline  jet  as  was 
observed  in  Fig.  3.2.20.  In  the  Mj=1.4  baseline  jet,  the  Strouhal  number  at  the  measured  fundamental 
sereeeh  frequeney  (fs)  is  about  0.37  (=  fsD/Uj)  as  shown  in  Fig.  3.2.2.  This  sereeeh  and  the  periodie 
struetures  seen  in  Fig.  3.2.20f  were  generated  by  the  feedbaek  meehanism.  As  in  the  over-expanded  jet. 
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the  initial  shear  layer  at  Mj  =  1 .4  is  exposed  to  two  sourees  of  perturbation:  one  is  seeded  by  the  plasma 
aetuators  and  the  other  is  amplified  by  the  flow-aeoustie  feedbaek  meehanism.  Thus,  Figures  3.2.20b  and 
3.2.21  suggest  that  there  is  a  strong  eompetition  for  energy  between  these  two  sets  of  struetures, 
espeeially  at  Stop  numbers  near  0.3  and  1.0.  When  the  foreed  struetures  are  not  energetie  either  due  to 
laek  of  organization  (at  low  Stop  numbers)  or  small  in  dimension  (at  high  Stop  numbers),  the  naturally 
amplified  struetures  prevail  over  foreed  struetures.  It  seems  that  the  jet  responds  either  to  the  foreing  by 
the  aetuators  or  to  the  natural  perturbation  by  the  flow-aeoustie  feedbaek  meehanism.  The  strueture 
spaeing  is  different  from  that  of  the  baseline  if  the  jet  responds  to  the  perturbation  seeded  by  the  aetuators. 

The  jet  width  for  the  three  baseline/unforeed  jets  is  eompared  as  shown  in  Fig.  3.2.22.  The  jet 
width  is  inereased  by  the  naturally  amplified  struetures  in  the  over-  or  under-expanded  baseline  jets.  Thus, 
for  imperfeetly  expanded  jets,  the  redueed  mixing  enhaneement  over  eaeh  eorresponding  baseline  jet  is 
partially  related  to  the  already  enhaneed  mixing  eaused  by  the  feedbaek  meehanism  (shown  in  Figs.  3.2.8, 
3.2.10,  and  3.2.12).  At  high  Stop^s,  the  jet  growth  is  redueed  by  foreing  at  some  azimuthal  modes  (not 
shown).  As  diseussed  earlier,  the  reduetion  in  jet  growth  is  due  to  suppression  of  naturally  amplified 
struetures  for  the  imperfeetly-expanded  jets. 


Fig.  3.2.22  Jet  width  for  the  three  baseline  jets. 
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3.3  Mach  0.9  Heated  Subsonic  Jet 

The  performance  of  plasma  actuators  in  an  Mj  0.9  subsonic  jet  was  investigated  for  stagnation 
temperature  ratio  {TJTa)  of  1.0,  1.4,  and  2.0  at  a  fixed  forcing  Strouhal  number  {SIdf  =  0.3),  which  is 
close  to  the  jet  preferred  instability.  The  focus  in  this  section  is  determining  whether  the  performance  of 
the  plasma  actuators  in  the  heated  jet  is  comparable  to  the  corresponding  cold  jet  for  limited  azimuthal 
modes  of  m  =  0,  1,  and  ±1.  In  this  section,  FWHM  indicates  the  jet  width  while  5  was  used  in  earlier 
sections. 

3.3.1  Mean  Flow  Results 

The  mean  flow  results  use  centerline  Mach  number  and  turbulent  kinetic  energy  as  well  as  jet 
width  as  the  metrics  for  mixing  enhancement.  It  should  be  noted  that  the  plotting  markers,  created  by 
severely  down  sampling  the  data,  are  for  visual  identification. 


Fig.  3.3.1  Centerline  Mach  number  for  jet  with  various  temperatures  and  forcing  azimuthal  modes. 


Figure  3.3.1  shows  the  centerline  Mach  number  for  all  the  cases  experimented.  Mach  number  for 
the  three  temperature  ratios  is  offset  in  increments  of  0.3  to  facilitate  visual  examination.  These  results 
show  that  the  jet  potential  core  ends  (the  centerline  Mach  number  starts  dropping)  at  smaller  x/D  across 
all  cases  as  temperature  is  increased.  In  all  forcing  cases,  the  jet  potential  core  is  shortened  by  forcing, 
especially  for  m  =  0  at  higher  temperatures,  with  mode  ±1  consistently  producing  the  greatest  centerline 
Mach  number  decay  rate.  It  should  be  noted  that  the  forcing  effectiveness,  indicated  in  this  figure  by  the 
rate  of  centerline  decay,  appears  to  be  improving  as  the  temperature  is  raised.  There  exists  a  small,  rapid, 
decay  region  upstream  of  x/D  =  4  for  m  =  0  which  becomes  stronger  as  the  temperature  increases.  It  is 
barely  observable  at  TJTa  =  1.4,  but  is  readily  apparent  at  2.0.  This  feature  most  likely  is  due  to  two 
events.  First,  Suzuki  and  Colonius  [2006]  have  shown  in  a  recent  work  that  m  =  0  mode  growth  in  an 
unforced  jet  is  significantly  more  enhanced  by  heating  of  the  jet  than  that  of  modes  1  and  2.  Second,  the 
nature  of  the  generated  vortex  pattern  and  interaction  of  coherent  vortices  in  the  axisymmetric  m  =  0  case 
creates  self-induction  which  generates  regions  of  additionally  accelerated  flow  as  can  be  clearly  seen  in 
To/Ta=  2.0  on  a  streamwise  range  of  3  to  4  x/D  [Kim  et  al.  2009a].  Additionally,  the  effectiveness  of 
mode  0  excitation  and,  to  a  lesser  extent,  mode  1  increase  significantly  in  the  elevated  temperature  cases 
and  their  downstream  centerline  Mach  number  profiles  (x/D  >  8)  become  similar.  These  results  indicate 
that  mixing  enhancement  due  to  forcing  is  achieved  at  all  temperatures  and  that  mixing  enhancement 
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improves  as  temperature  inereases.  Additionally,  the  results  indieate  a  signifieant  ehange  in  jet 
eharaeteristies  as  a  fiinetion  of  temperature  supported  by  the  signifieant  ehanges  in  m  =  0  mode  response. 


Fig.  3.3.2  Ratio  of  centerline  Mach  number  for  forced  jet  to  the  baseline  centerline  Mach  number. 


In  Fig.  3.3.2,  the  eenterline  Maeh  number  for  eaeh  mode  and  temperature  ratio  has  been  divided 
by  the  baseline  eenterline  Maeh  number  at  that  temperature.  This  provides  a  elear  view  of  the  profile 
ehanges  due  to  exeitation  (a  value  of  less  than  one  indieates  an  inerease  in  eenterline  Maeh  number 
deeay).  In  all  three  foreing  eases,  a  eonsistent  inerease  in  eenterline  deeay  is  observed.  The  eonvergenee 
by  x/D  =  12  of  the  two  elevated  temperatures  indieates  the  end  of  heating  effeets,  due  to  enhaneed  mixing 
and  signifieant  entrainment  of  the  ambient  air  into  the  jet.  This  figure  eliminates  the  ambiguity  in  the 
previous  figure  showing  definite  improvement  in  eontrol  effeetiveness  aeross  all  modes  as  temperature  is 
inereased.  It  also  elearly  shows  that  the  m  =  ±1  mode  is  produeing  the  greatest  inerease  in  eenterline 
deeay  rate. 

The  jet  width  eharaeteristies  are  shown  in  Fig.  3.3.3.  From  Fig.  3.3.3a,  there  are  a  few  notieeable 
ehanges  in  the  unforeed  spreading  of  the  jet  with  temperature.  As  the  temperature  is  raised,  the  spreading 
rate,  whieh  is  the  slope  of  the  eurves  shown,  slows  down  signifieantly  in  early  jet  development,  but 
beeomes  similar  beyond  x/D  ~  8.  On  average,  the  jet  with  temperature  ratio  of  2.0  is  about  9%  narrower 
than  the  unheated  jet.  At  a  temperature  ratio  of  2.0,  the  jet  has  no  signifieant  spreading  until  about  x/D  =  5. 
In  all  three  temperatures,  there  are  two  distinet  spreading  rates.  The  transition  between  the  two  spreading 
rates  oeeurs  between  x/D  =  5  and  7.  While  the  initial  spreading  rates  vary  substantially  (-0.03FWHM/x 
with  variation  of  45%),  the  seeond  spreading  rate  {-O.llFWHM/x  with  variation  of  3.5%)  is  about  the 
same  for  all  three  temperatures. 

It  may  be  eonfusing  to  some  readers  how  the  jet  potential  eore  ean  beeome  shorter  while  the  jet  is 
also  narrowing.  There  are  several  faetors  whieh  eontribute  to  this  result.  The  shortening  of  the  potential 
eore  due  to  heating  has  also  been  reported  in  Lepieovsky’s  work  [1999].  In  this  paper,  Lepieovsky  points 
out  that  as  temperature  inereases,  kinematie  viseosity  and  jet  exit  veloeity  inerease.  The  ehange  in 
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viscosity  affects  the  growth  of  the  shear  layer.  The  increasing  disparity  between  the  core  and  ambient 
velocities  increases  the  instability  of  the  shear  layer  which  gives  rise  to  vortices.  Secondly,  it  is  well 
known,  but  often  overlooked,  that  shear  layers  do  not  spread  symmetrically  about  the  plane  of  their  origin. 
The  exact  nature  of  the  spreading  dynamics  is  beyond  the  scope  of  this  paper,  but  the  work  of  Elliott  and 
Samimy  [1990]  is  a  representative  example  of  the  understanding  in  this  area.  Given  this  asymmetric 
behavior,  it  is  reasonable  to  accept  that  the  potential  core  can  shorten  while  the  jet  narrows. 


a)  Jet  width  in  an  unforced  jet  for  various  b)  A(Jet  spreading)  for  various  forced  cases, 

temperatures. 

Fig.  3.3.3  Jet  full  width  at  half  maximum  for  various  conditions. 


In  Fig.  3.3.3b  the  appropriate  baseline  jet  width  is  subtracted  from  the  forced  case  -  the  results  for 
temperature  ratios  of  1.4  and  2.0  are  artificially  shifted  up  for  visual  clarity.  It  should  be  noted  that, 
because  the  m  =  ±1  mode  is  elliptical  and  the  measurement  plane  is  on  the  major  axis  of  the  elliptic  jet 
cross-section,  the  data  collected  for  this  mode  exaggerates  the  increase  in  spreading  compared  to  the  total 
area  of  the  jet  [Kim  et  al.  2009a  or  Sec.  3.1].  In  the  previous  work,  it  is  noted  that  the  jet  width  on  the 
non-flapping  axis  can  be  approximated  by  the  unforced  width.  Consequently  the  width  data  for  the  m  = 
±1  mode  is  corrected  to  reflect  the  area  of  an  axisymmetric  jet.  This  correction,  known  as  the  equivalent 

jet  width  {deq),  is  calculated  as  5^^  =  where  ‘‘Base''  and  “±\"  are  the  measured  unforced  and 

forced  widths  respectively.  The  other  modes  are  axisymmetric  so  this  consideration  is  not  necessary.  In 
the  unheated  jet,  modes  0  and  1  have  distinct  spreading  rates,  but  in  the  elevated  temperature  cases,  their 
spreading  patterns  are  very  similar.  For  modes  0  and  1,  additional  spreading  due  to  forcing  is  not 
observed  until  about  x/D  =  5  or  6  indicating  that  forcing  at  these  conditions  does  not  have  an  impact  until 
the  jet  undergoes  the  spreading  transition  observed  in  the  unforced  cases.  Furthermore,  the  additional 
spreading  generated  by  forcing  with  modes  0  and  1  is  fairly  minimal.  On  the  other  hand,  mode  ±1  creates 
substantial  additional  spreading  with  onset  occurring  at  a  much  lower  x/D{-  3). 
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(a)  7/7,  =  1.0 


(b)  7/7,  =  1.4 


(c)  7/7,  =  2.0 

Fig.  3.3.4  Jet  centerline  turbulent  kinetic  energy  for  various  temperature  ratios  and  forcing 
modes. 


The  normalized  TKE  results  shown  in  Fig.  3.3.4  reinforee  the  results  observed  above.  It  should  be 
noted  that  TKE  is  a  two  eomponent  kinetie  energy  as  only  streamwise  and  normal  to  streamwise  veloeity 
eomponents  were  measured.  All  foreing  eases  show  signifieantly  higher  TKE  growth  and  they  saturate 
before  the  baseline  ease  does,  with  m  =  ±\  eonsistently  showing  the  greatest  peak.  The  TJTa=  1.0  ease  is 
eonsistent  with  results  obtained  previously  at  GDTL  [Samimy  et  al.  2007b].  As  the  temperature  is  raised, 
the  baseline  remains  relatively  unehanged  while  the  foreing  eases  grow  and  saturate  more  rapidly.  There 
is  little  ehange  in  the  m  =  0  and  m  =  1  modes  between  1.4  and  2.0,  but  m  =  ±1  ease  has  a  notieeable  peak 
inerease.  The  TKE  growth  onset  beeomes  more  rapid  in  all  foreing  eases  as  temperature  inereases.  The  m 
=  0  mode  shows  two  saturation  events  at  elevated  temperatures  supporting  the  assertion  of  the  presenee  of 
strong,  rapidly  developing,  struetures.  While  the  baseline  onset  does  also  beeome  more  rapid,  it  is  not 
nearly  as  dramatie  as  the  foreing  eases.  These  results  further  support  the  eonelusion  that  signifieant 
mixing  enhaneement  oeeurs  at  all  temperatures  with  inereased  mixing  enhaneement  at  higher 
temperatures. 
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3.3.2  Conditionally-Averaged  Results 

Now  that  the  mean  flow  eharaeteristies  have  been  examined,  the  stmetures  present  in  the  flow 
will  be  diseussed  in  an  effort  to  better  understand  the  effeets  of  foreing  on  the  flow  strueture  in  a  heated 
jet  and  explain  the  underlying  behavior. 

Conditional-averaging  is  eondueted  using  a  eorrelation  map  to  reveal  dominant  feature  of  the 
flow.  The  eonditional-averaging  in  this  seetion  is  different  from  that  in  earlier  seetions,  but  largely  very 
similar.  Eaeh  image  is  eorrelated  with  all  other  images.  If  the  stmetures  in  the  flow  are  relatively  eoherent, 
there  will  be  signifieant  eorrelation  among  a  subset  of  the  images  eolleeted.  Onee  the  eonditionally- 
averaged  veloeity  field  is  generated,  various  flow  analysis  methodologies  ean  be  applied  to  examine  the 
stmetures. 

In  the  data  presented  in  this  paper,  the  fiuetuating  v-eomponent  tended  to  bring  out  the  most 
readily  identifiable  stmetures  so  it  is  used  for  the  eorrelations.  The  region  of  interest  is  the  mixing  layer 
surrounding  the  jet  potential  eore  -  thus  eorrelation  is  performed  on  a  region  of  typieal  size  -1  to  +1  y/D 
by  0  to  4  x/D.  A  requirement  of  10%  usage  of  the  data  set  (typieally  about  70  fields)  for  eonditional- 
averaging  ensures  that  the  results  are  eonverged  and  the  stmetures  analyzed  are  statistieally  signifieant.  It 
should  be  noted  that  eonditional-averaging  of  the  Baseline  (no  foreing)  veloeity  fields  is  eonsiderably  less 
reliable  than  the  foreed  eases,  as  the  stmetures  in  the  flow  are  not  suffieiently  eoherent.  While  the  foreed 
eases  tend  to  satisfy  the  10%  requirement  with  a  eorrelation  eoeffieient  of  -0.55,  the  Baseline  eases 
typieally  only  meet  the  requirement  at  a  value  of  about  0.22.  This  is  expeeted  sinee  the  foreing  generates 
large  repeatable  stmetures  whereas  the  Baseline  jet  possesses  large  seale  but  not  quite  eoherent  stmetures. 

Figures  3.3.5  and  3.3.6  show  the  eonditionally-averaged  streamwise  eomponent  of  veloeity. 
These  figures  are  presented  to  give  a  visual  sense  for  the  flow  before  diseussing  the  stmetures.  In  the 
Baseline,  small,  transverse  oseillations  in  the  potential  eore  are  just  barely  observable.  The  small 
magnitude  of  these  oseillations  indieates  a  laek  of  eoherenee  in  large-seale  stmetures.  The  results  for  eaeh 
foreed  ease  show  signifieant  alteration  of  the  jet  potential  eore.  In  the  m  =  1  mode,  there  are  pronouneed, 
transverse  oseillations  and  the  potential  eore  is  shortened  by  about  2  x/D  eompared  to  the  baseline.  The  m 
=  ±1  mode  exhibits  qualitatively  similar  results  to  mode  1,  but  the  oseillation  amplitude  is  greater  and  the 
potential  eore  is  shortened  by  3  x/D.  In  the  m  =  ±1  mode,  the  oseillations  are  so  large  that  (as  seen  in  Fig. 
3.3.5d  at  x/D  =  3)  the  potential  eore  is  perturbed  entirely  to  one  side  of  the  jet  eenterline.  As  was 
diseussed  in  the  mean  flow  results,  the  baseline,  m  =  1,  and  m  =  ±1  modes  vary  only  quantitatively  as  the 
temperature  is  raised  -  the  potential  eore  beeomes  shorter. 

The  behavior  of  the  m  =  0  mode  is  more  eomplex  than  the  others  as  ean  be  seen  in  Figs.  3.3.5b 
and  3.3.6.  In  the  unheated  ease,  there  are  weak  longitudinal  variations  in  the  potential  eore  whieh  are 
aeeompanied  by  axisymmetrie  variation  (waists)  whieh  is  out  of  phase  relative  to  the  longitudinal 
variations.  Referring  to  the  mean  flow  results,  there  is  almost  no  ehange  in  the  length  of  potential  eore  due 
to  foreing  at  this  mode  and  temperature  (Fig.  3.3.1).  As  the  temperature  inereases,  the  amplitudes  of  both 
the  transverse  and  longitudinal  variations  inerease.  The  low  veloeity  waist  regions  beeome  signifieantly 
more  pronouneed.  Additionally,  the  effeet  of  foreing  with  this  mode  on  the  potential  eore  length  inereases 
substantially.  Combining  inspeetion  of  Figs.  3.3.1,  3.3.5b,  and  3.3.6,  the  deerease  in  potential  eore  length 
due  to  m  =  0  foreing  is  about:  1,  4,  and  4.5  x/D  in  order  of  inereasing  temperature.  These  results  support  a 
signifieant  ehange  in  jet  eharaeteristies  as  was  diseussed  in  the  mean  flow  results  (Fig.  3.3.2). 
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Fig.  3.3.5  Conditionally  averaged  streamwise  velocity  fields  for  T/Ta  =  1.0. 
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Fig.  3.3.6  Conditionally  averaged  streamwise  velocity  fields  for  hi  =  0. 
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Fig.  3.3.7  Swirling  strength  (Hz)  and  Galilean  streamlines  for  T/T^  =  1.0. 


Figures  3.3.7  and  3.3.8  show  the  eonditionally-averaged  images  using  two  vortex  identifieation 
methods.  The  primary  tool  used  to  diseuss  struetures  in  this  paper  is  known  as  the  swirling  strength 
method  -  originally  developed  by  Chong  et  al.  [1990].  This  vortex  identifieation  method  uses  the 
imaginary  eomponent  of  the  eigenvalues  of  the  veloeity  gradient  tensor  to  loeate  regions  of  rotating  fluid 
while  ignoring  regions  of  pure  shear.  Compared  to  vortieity,  swirling  strength  is  a  superior  tool  in  flow 
strueture  analysis  due  to  the  elimination  of  sheared  flow  eontributions.  For  two-dimensional  imaging,  the 
essenee  of  this  teehnique  is  that  regions  of  swirl  obey  the  following  eondition  for  the  produet  of  the  two 
shear  eomponents. 


< 


(3.3.1) 


where  u  and  v  are  the  x-  and  y-eomponents  of  the  veloeity,  respeetively,  and  a  subseript  indieates  a 
derivative  with  respeet  to  that  variable.  The  most  important  result  of  this  inequality  is  that  the  imaginary 
eomponent  of  the  eigenvalues  will  be  non-zero  only  if  this  produet  is  negative,  whieh  only  oeeurs  in 
regions  where  the  flow  is  rotating.  The  imaginary  eomponent  is  shown  as  a  eontour  plot  in  the  following 
figures  and  has  units  of  inverse  seeonds  -  swirling  strength  is  a  measure  of  the  rotational  rate  of  the  fluid 
strueture. 
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The  peaks  in  swirling  strength  are  used  to  ealeulate  several  empirieal  values  for  deseribing  the 
struetures  in  the  flow  -  these  values  are  shown  in  Table  4.  The  peak  values  for  the  set  of  struetures  in  a 
given  plot  are  averaged  to  produee  the  Peak  Swirling  Strength  (PSS).  The  U-eomponent  of  the  veloeity  at 
eaeh  peak  is  averaged  to  determine  the  average  eonveetive  veloeity  (f/c).  This  eonveetive  veloeity  is  used 
to  determine  the  Galilean  deeomposition  and  to  produee  the  streamlines  shown  in  the  figures.  Lastly,  the 
peak  eoordinates  are  used  to  ealeulate  the  spaeing  between  struetures. 


a)  777,  =  1.4  b)  777,  =  2.0 

Fig.  3.3.8  Swirling  strength  (Hz)  and  Galilean  streamlines  for  /w  =  0  in  heated  jet. 


The  flow  strueture  plots  for  the  unheated  jet  are  shown  in  Fig.  3.3.7.  The  Galilean  streamlines 
show  good  agreement  with  the  identified  swirling  regions,  and  also  provide  additional  information  on  the 
shape  of  the  struetures.  Furthermore,  the  streamlines  indieate  that  some  swirling  regions  are  not  vortex 
eenters.  These  regions  are  most  likely  non-physieal  artifaets  of  finite  PIV  grid  resolution  and  the  low 
number  of  images  used  in  eonditional-averaging.  As  expeeted,  the  struetures  in  the  foreed  eases  are  more 
organized,  and  are  apparent  almost  immediately  downstream  of  nozzle  exit.  Baseline  struetures  are  less 
organized,  flattened,  and  don’t  beeome  signifieantly  organized  until  farther  downstream  eompared  to 
those  in  the  foreed  eases.  The  symmetrie  nature  of  mode  0  is  immediately  apparent.  The  other  two  modes 
have  a  similar  appearanee,  but  their  distinetions  are  borne  out  by  the  data  in  Table  4.  The  similarity  of  the 
foreed  ease  strueture  spaeing  to  the  unforeed  eases  shows  that  the  frequeney  ehosen  is  indeed  elose  to  the 
natural  instability  frequeney  of  the  jet.  All  of  the  foreed  eases  have  higher  PSS  than  the  baseline.  Mode  ±1 
is  observed  to  have  the  strongest  PSS. 

The  similarity  of  the  baseline  and  m  =  1  strueture  patterns  warrants  additional  diseussion.  First,  as 
mentioned  above,  the  eorrelation  level  of  the  baseline  veloeity  field  (-  0.22)  is  mueh  lower  than  foreed 
eases  (-0.55)  so  the  struetures  analyzed  are  less  representative  of  the  average  jet  behavior.  Next,  previous 
work  by  Tinney  et  al.  [2008]  on  the  azimuthal  mode  eomposition  of  an  unforeed  Maeh  0.85  jet  shows  that 
low  order  helieal  modes  are  dominant  in  the  jet.  All  azimuthal  modes  produee  a  strueture  pattern  whieh  is 
either  staggered  or  symmetrie  -  any  single  helieal  mode  will  produee  a  staggered  strueture  pattern.  Sinee 
eonditional-averaging  will  seek  out  a  pattern  no  matter  how  weak,  it  is  reasonable  to  expeet  the  algorithm 
to  find  one  of  these  patterns.  Additionally,  in  a  study  eondueted  by  Moore  [1977]  in  high  subsonie  jets. 
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large  scale  structures  in  an  unforced  jet  were  observed.  By  synchronizing  schlieren  images  with  the 
arrival  of  structures  on  one  side  of  the  jet  as  determined  by  a  microphone,  he  was  able  to  create  a 
conditionally  averaged  image  of  the  jet  structures.  What  he  found  was  that  the  structures  on  the  other  side 
of  the  jet  were  absent  from  this  averaged  image.  When  he  examined  individual  frames  he  saw  that  the 
structures  were  sometimes  symmetrically  arranged  and  in  other  frames  they  were  staggered.  Since  the 
conditional-averaging  used  in  the  present  paper  operates  over  the  whole  jet  cross-section,  it  makes  sense 
that  one  structural  pattern  would  be  isolated.  Lastly,  while  the  baseline  structures  shown  in  Fig.  3.3.7  may 
not  be  representative  of  the  jet  in  an  average  sense,  they  are  (by  the  nature  of  the  conditional-averaging 
routine)  representative  of  the  strongest  structures  present  in  the  jet. 

As  before,  the  mode  0  case  is  presented  alone  in  Fig.  3.3.8  due  to  the  general  similarity  of  other 
cases  over  the  temperature  range.  Comparing  Fig.  3.3.7b  and  Fig.  3.3.8,  it  is  apparent  that,  as  the 
temperature  increases,  the  organized  structures  become  confined  to  a  smaller  region  closer  to  the  nozzle 
exit.  The  decreasing  region  of  organized  structures  reinforces  the  conclusions  drawn  from  the  mean  flow 
and  velocity  results  discussed  in  the  previous  sections  -  a  significant  change  takes  place  in  jet 
characteristics  with  increasing  temperature. 

When  the  temperature  ratio  is  raised  to  1.4,  there  are  several  observable  changes  in  the  flow 
structures.  The  basic  flow  patterns  of  each  of  the  four  cases  stay  the  same,  but  the  amplitude,  spacing,  and 
streamwise  extent  of  the  structures  change.  The  PSS  increases  for  each  case.  The  PSS  of  the  forced  cases 
continues  to  be  significantly  higher  than  the  baseline  case.  The  most  interesting  change  is  that  the  PSS  of 
mode  0  is  now  slightly  greater  than  that  of  mode  ±1.  At  the  temperature  ratio  of  2.0,  the  same  trends 
observed  in  1.4  continue  and  now  mode  0  clearly  has  the  largest  PSS. 

The  behavior  of  the  convective  velocity  merits  some  discussion.  According  to  theoretical 
equations  for  convective  velocity  such  as  those  cited  in  Thurow  et  al.  [2008],  the  non-dimensionalized 
convective  velocities  for  the  three  temperature  ratios  under  discussion  are:  0.53,  0.48,  and  0.44  (in 
ascending  order  of  temperature).  While  there  is  fairly  good  agreement  with  theory  in  the  unheated  cases, 
there  is  a  significant  departure  from  theory  in  the  elevated  cases.  According  to  theory,  the  ratio  of 
convective  velocity  to  exit  velocity  should  decrease  with  increasing  temperature,  but  the  data  collected 
shows  that  this  ratio  is  increasing,  except  for  the  m  =  0  case.  Soteriou  and  Ghoniem  [1995]  observed  that 
convective  velocity  decreased  as  density  ratio  decreased.  In  those  experiments,  the  velocity  ratio  (jet 
velocity  to  ambient  velocity)  was  constant  (-1.5).  Performing  a  linear  fit  with  temperature  ratio  as  the 
independent  variable,  the  slope  of  the  velocity  ratio,  using  a  small  constant  velocity  for  the  ambient  flow, 
in  these  experiments  (-23)  is  greater  than  the  slope  of  the  ratio  of  convective  velocity  to  the  ambient 
(~15).  It  is  concluded  that  the  convective  velocity  is,  qualitatively,  in  agreement  with  Soteriou  and 
Ghoniem’ s  observations.  The  decrease  in  non-dimensionalized  convective  velocity  in  mode  0  at  a 
temperature  ratio  of  1 .4  cannot  be  explained  at  this  time,  but  a  likely  suspect  is  the  significant  increase  in 
structure  strength.  Unfortunately,  the  extreme  environment  of  heated,  high  Reynolds  number  jets  makes  it 
difficult  to  acquire  data  on  this  subject,  which  limits  the  number  of  papers  published  in  the  literature  on 
the  subject  and,  consequently,  the  conclusions  which  can  be  drawn  here. 

The  trends  in  PSS  are  made  clearer  by  non-dimensionalizing  the  data.  Generally,  the  strength  of 
the  vortices  is  directly  correlated  to  the  jet  exit  velocity  as  evidenced  by  the  consistency  of  the  non- 
dimensional  PSS.  This  indicates  that  the  forcing  mechanism  on  average  is  not  changing  as  the 
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temperature  ehanges  -  at  least  for  this  Strouhal  number.  The  notable  exeeption  to  this  eonsisteney  is  mode 
0.  At  a  temperature  ratio  of  1.0,  mode  0  strength  is  eomparable  to  mode  1.  At  the  temperature  ratio  1.4, 
the  strength  of  mode  0  has  risen  and  is  eomparable  to  mode  ±1.  This  suggests  a  ehange  in  the  reeeptivity 
of  the  jet  to  foreing  at  this  mode  and  Strouhal  number.  As  was  diseussed  previously,  inereasing  mode  0 
energy  with  inereasing  temperature  relative  to  other  modes  has  been  observed  [Suzuki  and  Colonius  2006] 
so  it  is  reasonable  to  expeet  an  inerease  in  strueture  strength. 


Table  4  Structure  information  for  all  cases. 


Mode 

TJT„ 

ud  Uj 

Conditionally 

Peak  Swirling 

D*PSS 

Averaged:  A/D 

Strength  (Hz) 

Uj 

Baseline 

1.0 

0.57 

1.61 

119 

0.011 

Baseline 

1.4 

0.58 

1.55 

130 

0.010 

Baseline 

2.0 

0.55 

1.36 

168 

0.011 

0 

1.0 

0.51 

1.55 

177 

0.016 

0 

1.4 

0.45 

1.46 

286 

0.022 

0 

2.0 

0.54 

1.36 

324 

0.020 

1 

1.0 

0.52 

1.67 

175 

0.016 

1 

1.4 

0.52 

1.52 

221 

0.017 

1 

2.0 

0.60 

1.46 

249 

0.016 

±1 

1.0 

0.54 

1.65 

221 

0.020 

±1 

1.4 

0.56 

1.62 

270 

0.020 

±1 

2.0 

0.59 

1.49 

294 

0.019 

A  decrease  in  structure  spacing  is  expected  because  of  the  changing  ratio  of  convective  velocity 

to  exit  velocity/^  =  ^ j j^^df  )  •  Using  the  theoretical  values  for  convective  velocity,  the  structure 

spacing  is  1.66,  1.50,  and  1.36  x/D  in  order  of  increasing  temperature  ratio.  While  at  first  this  may  look 
promisingly  related  to  the  results  in  Table  4,  the  data  resolution  must  be  remembered.  The  range  of 
variation  in  the  data  is  0.31  x/D,  but  the  data  resolution  is  0.11  x/D.  Due  to  low  resolution  relative  to  the 
range  of  variation  in  structure  spacing,  the  only  valid  conclusion  is  that  the  structure  spacing  is 
approximately  the  same  for  all  cases.  The  results  for  the  Mj  1.3  unheated  jet  confirms  that  the  structure 
spacing  is  almost  the  same  for  all  azimuthal  modes  when  the  forcing  frequency  is  the  same  (Fig.  3.2.16). 
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3.4  Mach  1.65  Cold  Jet 

Limited  experiments  are  eondueted  for  Mj  1.65  unheated  jets  from  eontoured  and  eonieal  nozzles 
(referred  as  eontoured  and  eonieal  jet,  respeetively).  The  foeus  of  this  seetion  is  to  find  whether  the 
LAPP  As  have  eontrol  authority  in  a  supersonie  jet  with  a  high  jet  Maeh  number  (Mj  =  1.65)  and  in  a 
shoek  eontaining  eonieal  jet.  Based  on  the  results  in  Mj  0.9  and  1.3  jets,  azimuthal  mode  and  SIdf’s  are 
limited  to  m  =  ±1  and  Stop  =  0.18,  0.29,  0.47,  and  1.4.  Note  that  the  results  for  Mj  1.65  jets  are  very 
preliminary  and  investigation  is  still  underway. 


Shown  in  Fig.  3.4.1  are  Maeh  number  and  TKE  along  the  eenterline  of  Mj  1.65  eontoured  and 
eonieal  jets.  The  eenterline  Maeh  number  for  the  eontoured  jet  shows  slight  undulation  although  it  was 
operated  in  the  design  eondition  as  in  the  Mj  1.3  perfeetly  expanded  jet  (See.  3.2).  For  the  eonieal  jet  (Fig. 
3.4. le),  the  eenterline  Maeh  number  undulates  signifieantly  due  to  stronger  shoek  eell  trains  in  the  jet 
plume  eompared  to  the  eontoured  jet  as  manifested  in  sehlieren  images  in  Fig.  3.4.2.  The  strong  shoek 
trains  are  expeeted  sinee  the  eonieal  nozzle  has  a  sharp  throat  and  a  eonieal  diverging  seetion. 


x/D 

(a)  Centerline  Mach  number  in  the  contoured  jet 


x/D 

(c)  Centerline  Mach  number  in  the  conical  jet 


x/D 

(b)  TKE  along  the  jet  centerline  in  the  contoured  jet 


x/D 

(d)  TKE  along  the  jet  centerline  in  the  conical  jet 


Fig.  3.4.1  Mach  number  and  TKE  along  the  jet  centerline  in  m  =  ±1 


It  seems  that  the  eenterline  Maeh  number  deeay  in  the  Mj  1.65  eontoured  jet  is  qualitatively 
similar  to  that  in  Mj  0.9  subsonie  (Fig.  3.1.7)  and  Mj  1.3  perfeetly-expanded  (Fig.  3.2.8)  jets,  although 
eomparison  is  not  straightforward  sinee  the  potential  eore  length  and  Maeh  number  are  not  the  same.  The 
optimum  Stop,  produeing  maximum  jet  spreading,  is  about  0.3  similar  to  those  in  Mj  0.9  and  1.3  jets.  The 
effeets  of  foreing  deerease  as  the  Stop  inereases  -  also  similar  to  what  was  observed  in  Mj  0.9  and  1.3  jets. 
The  TKE  level  signifieantly  inereases  at  low  SIdf^s,  but  the  saturated  TKE  level  is  mueh  lower  than  those 
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observed  in  Mj  0.9  and  1.3  jets  (Figs.  3.1.9  and  3.2.13,  respeetively).  This  implies  that  the  foreing  in  Mj 
1.65  eontoured  jet  is  slightly  less  effeetive  than  that  in  low  Mj  jets. 

It  is  expeeted  that  foreing  in  the  Mj  1.65  eonieal  jet  to  be  similar  to  that  in  Mj  1.2  overexpanded 
and  Mj  1.4  underexpanded  jets  (See.  3.2),  sinee  there  will  be  naturally-amplified  struetures  in  the  Mj  1.65 
eonieal  jet  as  in  the  over-  and  under-expanded  jets.  In  Mj  1.2  overexpanded  and  Mj  1.4  underexpanded 
jets,  the  eenterline  Maeh  number  deeay  was  not  enhaneed  by  foreing  (Figs.  3.2.10  and  3.2.12).  However, 
the  eenterline  Maeh  number  deeay  is  signifieantly  enhaneed  in  the  Mj  1.65  eonieal  jet,  espeeially  at  Stop  ^ 
0.29  (Fig.  3.4. le).  As  a  result,  the  TKE  level  is  also  signifieantly  inereased  as  shown  in  Fig.  3.4. Id.  The 
effeets  of  foreing  at  other  Stores  are  limited  or  minimal  as  seen  in  eenterline  Maeh  number  and  TKE.  At  a 
high  Stop,  the  eenterline  Maeh  number  deeay  is  retarded  and  TKE  level  is  redueed  by  foreing.  Again,  this 
is  similar  to  what  was  observed  in  Mj  1.2  overexpanded  and  Mj  1.4  underexpanded  jets  (See.  3.2).  This 
similarity  suggests  that  the  foreed  struetures  prevail  at  intermediate  Stop’s  while  naturally-amplified 
struetures  do  at  low  and  high  Stop’s  as  ean  be  inferred  from  Figs.  3.2.19  and  3.2.21.  The  eompetition 
between  these  two  types  of  struetures  seems  to  be  responsible  for  the  redueed  effeetiveness  of  foreing  in 
the  eonieal  jet. 


(a)  Contoured  j  et  (b)  Conical  j  et 

3.4.2  Time-averaged  schlieren  images  for  the  baseline/unforced  jets. 
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4.  CONCLUDING  REMARKS 

Active  flow  control  of  jets  with  Localized  Arc  Filament  Plasma  Actuators  (LAPP As)  was 
conducted  over  a  wide  range  of  the  fully  expanded  jet  Mach  numbers  (Mj).  The  jet  Mach  numbers 
covered  in  the  present  research  are  0.9  (with  a  converging  nozzle),  1.2  (overexpanded),  1.3  (perfectly 
expanded),  and  1.4  (underexpanded)  with  a  design  Mach  number  1.3.  Additionally,  limited  experiments 
were  carried  out  for  an  Mj  =  1.65  perfectly-expanded  jet  with  both  contoured  and  conical  nozzles.  The 
exit  diameter  is  2.54  cm  (1  inch)  for  all  cases  and  eight  LAPP  As  are  equally  distributed  on  the  perimeter 
of  a  boron  nitride  nozzle  extension.  The  performance  of  LAPP  As  was  evaluated  using  the  centerline 
Mach  number  decay,  jet  width  and  TKE  developments,  calculated  from  two-dimensional  PIV  results. 

The  jet  spreading  was  strongly  dependent  on  duty  cycle,  forcing  frequency,  and  azimuthal  modes. 
The  effects  of  duty  cycle  was  investigated  at  m  =  ±1  for  several  forcing  frequencies  in  Mj  0.9  jet.  At  each 
fixed  forcing  frequency,  the  maximum  jet  spreading  occurred  at  the  smallest  duty  cycle  producing  stable 
plasma  arc  (referred  as  optimum  duty  cycle).  A  relationship  between  the  optimum  duty  cycle  and  forcing 
frequency  was  determined  from  the  extensive  results  in  the  Mj  0.9,  and  this  relation  was  used  for  all 
experiments.  In  Mj  0.9  subsonic  and  Mj  1.3  perfectly-expanded  jets,  very  extensive  experiments  were 
conducted  to  find  the  effects  of  forcing  frequency  and  azimuthal  modes.  The  forcing  Strouhal  number 
ranged  from  0.09  to  3.0  and  the  azimuthal  modes  investigated  were  m  =  0  -  3,  ±1,  ±2,  and  ±4,  only 
available  modes  with  eight  actuators.  The  performance  of  LAPP  As  was  also  dependent  on  the  stagnation 
temperature  of  the  jet  and  jet  Mach  number.  The  effects  of  stagnation  temperature  were  investigated  for 
1.0,  1.4,  and  2.0  the  ambient  temperature  only  in  Mj  0.9  jet  for  very  limited  azimuthal  modes  and  SIdf-  In 
an  Mj  1.65  perfectly-expanded  jet,  the  control  authority  of  LAPP  As  was  investigated  for  only  m  =  ±1  and 
Stop  0.3. 

The  jet  spreading  increased  with  decreasing  duty  cycle  as  long  as  there  were  a  minimal  number  of 
misfires  in  the  plasma  actuators.  Thus,  the  optimum  duty  cycle  for  each  forcing  frequency  was 
determined  to  be  the  lowest  value  which  ensured  stable  plasma  generation.  Extensive  experiments  in  Mj 
0.9  subsonic  and  1.3  perfectly-expanded  supersonic  jets  showed  that  the  jet  spreading  was  significantly 
enhanced  reaching  the  maximum  at  about  Stop  =  0.3  for  most  azimuthal  modes.  The  most  and  least 
effective  azimuthal  modes  were  m  =  ±1  and  ±4  respectively.  The  results  also  showed  that  the  effects  of 
forcing  were  very  similar  in  Mj  0.9  subsonic  and  1.3  perfectly-expanded  supersonic  jets  in  a  qualitative 
sense.  The  generated  structure  spacing  was  solely  determined  by  forcing  frequency,  while  the  strength 
and  pattern  of  structures  depended  on  both  forcing  frequency  and  azimuthal  mode. 

The  turbulent  kinetic  energy  along  the  jet  centerline  also  increased  significantly  near  Stop  =  0.3 
for  most  azimuthal  modes.  The  trend  of  TKE  amplification  was  similar  to  that  of  centerline  Mach  number 
decay,  implying  more  amplification  for  enhanced  jet  spreading.  Thus  the  maximum  TKE  amplification 
was  observed  in  m  =  ±1  and  at  a  Stop  near  0.3. 

The  results  in  the  heated  Mj  0.9  jet  showed  that  the  effects  of  forcing  increase  with  increasing 
stagnation  temperature.  In  addition,  the  jet  spreading  in  m  =  0  increased  relatively  more  significantly  and 
was  comparable  to  that  in  m  =  1  at  an  elevated  stagnation  temperature  while  it  was  one  of  the  less 
effective  modes  in  an  unheated  jet.  Possible  reasons  for  the  improved  actuation  effectiveness  are 
discussed 
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In  off-design  conditions  of  Mj  =1.2  (overexpanded)  and  1.4  (underexpanded),  the  forcing  was 
less  effective  compared  to  the  perfectly-expanded  case  of  Mj  =  1.3.  Flow  visualization,  via  Galilean 
streamlines,  shows  that  there  are  naturally-amplified  flow  structures  (generated  by  a  natural  feedback 
mechanism)  in  addition  to  the  structures  generated  by  forcing,  especially  at  low  and  high  forcing  Strouhal 
numbers.  The  competition  of  energy  between  these  two  structures  seemed  to  be  responsible  for  the 
reduced  effectiveness  of  forcing. 

The  performance  of  LAPP  As  in  a  high  Mj  supersonic  jet  (Mj  =  1.65)  was  investigated  at  several 
Stop’s  near  0.3  for  m  =  ±1.  The  results  showed  slightly  reduced  forcing  effectiveness  at  this  high  jet  Mach 
number  compared  to  those  in  Mj  0.9  and  1.3  jets.  It  has  not  been  determined  if  the  reduction  in 
effectiveness  is  due  to  a  lack  of  LAFPA  control  authority  or  increased  flow  compressibility. 
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